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SUMMARY 


The purpose of the work carried out under NASA Headquarters grant 
NGR 05-003-458 was two-fold: 

1) To examine models of Jupiter's magnetosphere with the particular 
purpose of predicting the X-ray flux that would be emitted in 
auroral or radiation zone processes. 

2) To investigate the various types of X-ray detection for energy 
resolution, efficiency, reliability, and background, and then 
from the model fluxes to decide under what models Jovian X-rays 
could be detected. 

We have found that the proportional counter is best suited for an outer 
planet mission because of its overall simplicity, the achievability of 
low background, the ease of attaining large areas at low mass with a 
counter of modest size, and its suitability for the MJS mission. It is 
concluded that there is a high probability that Jovian X-rays will be 
detected at distances up to 1 A.U. from the planet. 

The report consists of a discussion of the models used to predict 
X-ray fluxes and the instrumentation most suited to measure them, 
of the work reported here was carried out by Dr. Kevin Hurley. 


Much 
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FINAL REPORT FOR NGR 05-003-458 (SUPPORTING RESEARCH FOR 
OUTER PLANETS MISSION - EXPERIMENT DEFINITION PHASE) 

As the planet Jupiter is known to possess an active magnetosphere, 
and Saturn is strongly suspected to have one, these two planets should be 
sources of X-radiation, produced by bremsst rah lung from energetic electrons 
which precipitate into the planetary atmospheres. A study of planetary 
X-rays therefore gives information about particle fluxes, and constitutes 
a powerful tool for the remote sensing of magnet ospheric dynamics. We 
had made previous estimates of Jupiter's X-ray flux, based on theoretical 
considerations (1). Our goal, under NGR 05-003-458, was 

First, to revise the estimates of the Jovian X-ray fluxes, using 
any relevant new information, and to generate similar estimates for 
Saturn. 

Second, using these new estimates, to determine the best type of X-ray 
experiment, for a Grand Tour or MJS-type spacecraft, which could detect 
these fluxes. 

Third, to compare the best MJS experiment with other types of experi- 
ments, such as balloon, rocket, and satellite X-ray detectors, to determine 
whether a flyby experiment offered significant advantages over these other 
types . 

We have summarized our work on these three phases in the following 
three sections. - - ~ _ -- 
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I. MODELS FOR PLANETARY X-RAY EMISSION 

There is at present no data on the spectrum of kilovolt electrons 
in the Jovian magnetosphere. Therefore, in order to calculate the flux 
of Jovian X-rays, models must be constructed which allow the known terres- 
trial electron spectra to be scaled to the Jovian magnetosphere. Once 
the assumed Jovian electron spectrum has been calculated, the X-ray flux 
follows by applying the formulas for bremsstrahlung cross sections. This 
is explained in Appendix A. | 

MODEL IA 

i 

This model is based on the assumption that the ultimate source of 
energy for all terrestrial auroral phenomena is derived from the solar 
wind particle energy. If the same fraction of Jupiter's energy input 
goes into precipitating electron fluxes as the earth's, the Jovian 
fluxes can be found by scaling the terrestrial fluxes by a factor which 
depends on the intercepted solar wind energies. 

In the cases of both the earth and Jupiter, the planetary magnetic 
fields do not extend out into space indefinitely; they are balanced at 
some point by the pressure of the solar wind. Thus if n, m, and v are 
the number density, particle mass, and velocity of the solar wind, and 
B d is the planetary dipole field strength, the boundary will occur where 

, 2 2 ( 2 B d ) 2 

Znmv^cos^ y = — * — 

oTT 

Here, ¥ is the angle between v and a normal to the boundary, and the 
field just inside the boundary has been increased to twice the dipole 
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value, due to currents flowing along the boundary. Since 4 1 depends, in 
turn, upon the shape of the boundary, the above equation must be solved 
by iterative techniques, except at the subsolar point, where ¥ = 0. Thus 
if M is the planetary magnetic moment, and r is the distance from the 
center of the planet (assumed to be the center of dipole) , the radius 
of the magnetosphere at the subsolar point is 

r Q = (M 2 /4irnmv 2 ) V 6 . 


The quantity of interest is the cross sectional area presented to the 
solar wind, or the cross section in the dawn-dusk meridian. Detailed calcu 
lations (2) have shown that this cross section is approximately elliptical 
with semiminor and semimajor axes r Q and 1.4 r , respectively. Thus the 
cross sectional areas should scale as r Q 2 and the scaling factor relating 
Jovian to terrestrial intercepted energy becomes, in this case ; 


3 sn T mv 3 r 2 
J oJ 

Jjn mv d r 2 
t ofc 


M i R E 2 

m e r j 2 


2/3 


150 


where the subscripts J and E refer to Jupiter and the earth, and R >is 
the distance between the planet and the sun. It has been assumed that 

- I 

the solar wind velocity is the same at Jupiter as at the earth, while 
the number density decreases as 1/R 2 . Mj has been taken to be 4 x lO 30 
gauss-cm 3 , or~5 x io 4 M £ . r^ is. .then approximately 50 Rj. 

The total power intercepted by the earth's magnetosphere is about 
IO 20 ergs/sec, of which roughly— 10 18 - ergs/sec goes-into -auroral phenomena 
of all kinds. Jupiter, then, should intercept 1.5 * 10 22 ergs/sec and 
use 1.5 x 10 ergs/sec for its auroral phenomena, according to this 
model. (For comparison, the average power dissipated in Jovian 
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decimetric and decametric radiation is about 10 16 ergs/sec.) If the 
energy spectrum of precipitating electrons is now taken to be 

= C e E/E o electrons/cm 2 sec keV 
dfc 

(see Appendix A), the total energy flux precipitated is 

i 

E ^ dE = CE 2 keV/cm 2 sec . 

dE o - 

0 

A question remains as to the area into which electrons are precipitated. 
Although the earth's auroral zones are generally taken to be those areas 
between latitudes 65° and 70° in both the northern and southern hemispheres, 
(the boundary between open and closed field lines) , precipitation at any 
given time will involve only about 1% of this area, or 3 x 10 15 cm 2 over 
one hemisphere. If the same latitudes are assumed to define the Jovian 
auroral zones, and if the precipitation occurs over the same percentage 
of the area, the size of the precipitating region will scale as the square 
of the planetary radius. [However, see (3) for a different' approach 
which gives a slightly lower value.] Taking rj = 7.1 x 10 9 cm as the 
radius of Jupiter, and r^ = 6.3 x 10 8 cm as the radius of the earth, 
the precipitating region on Jupiter is found to be 110 times larger than 
the area on the earth, or 3 x 10 17 cm 2 . Although electron precipitation 
shows different time structure on the earth's day and night sides, the 
total X-ray energy is approximately’ the same for' both sides; thus these 
calculations should apply to both sides of a planet. The total power 
in electronifluxes becomes AjCjE q j 2 for Jupiter, and A^C^E^ 2 = 2.1 x 10 25 
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keV/sec for the earth, where Ag and Aj are the areas of the precipitating 
regions on the earth and Jupiter and the spectral parameters C and E q have 
been assumed to be different for the two planets. Now, scaling the power 
dissipated by the ratio of the intercepted energies, 

A t C t E 2 = 150 A_C„E T 2 = 3.15 x 10 27 keV/sec. 

J J oJ E E oJ 

I 

From Appendix A, Formula IA, the X-ray spectrum at the earth at MJS launch, 
when Jupiter is 8.25 x 10 13 cm away, is 


dh IA (ho)) 

d W 


= 8.83 x 10" 8 


-hco/Epj 

•hoi 


X-rays/cm 2 sec keV. 


This spectrum is plotted in Figure 1, for various values of the e -folding 

r- -I 

energy E q j. i 

MODEL IB 


One objection which can be raised to Model IA is that it is a steady 
state model: energy flows into and out of the magnetosphere without 

being stored. In reality, there is good evidence that the magnetic field 
lines in the earth's auroral zones connect- to a highly diamagnetic ;sheet 
of plasma which stretches out in the antisolar direction behind the earth 
(4, 5). Thus the electron energy density in this region is comparable to 
the magnetic field energy density, and energy is indeed stored in the tail 
of the magnetosphere. Assuming that this represents an upper limit for 
Jupiter as well as for the earth, the Jovian electron energy density may 
be scaled by the ratio of the magnetic field energy densities. If 

is the velocity of an electron with kinetic energy E and 



u 


c 



DNIR/DHU CPHDTCJNS/CM2 SEC KEU ) 
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and rest energy me 2 , the energy density of an exponential electron spectrum 

dN _ -E/E 0 ... ! 

— = C e ° is given by 

. oo 

f - C e" E/E ° dE = - E 3 / 2 r (3/2) 

J u c V 2 o 

0 

where r is the gamma function. The use of the nonrelativistic velocity 
is justified by the fact that, for the E o j to be used in this expression, 

the contribution at energies E ~ me 2 is very small. ! 

j 

j 

The scaling factor is Bj 2 /B^ 2 , where Bj and Bg are the magnetic 
field strengths at corresponding points on Jupiter and the earth; assuming 
dipole fields, this becomes I 




2 


2.5 x 10 3 . 


Thus 1 

C.E T 3 / 2 = 2.5 x 10 3 C_E 3 / 2 
J oJ E oE 

or 

C J E oJ 2a J ■ 2 - 5 * 103 c e e oE 3/2 a j ’ 1 - 66 " 1030 ■ 

* 

From Appendix A, formula. IA, the X-ray spectrum at the earth is 


dn IB (fioj) 

d(fuo)~ 


-fuo/E , 

4.66 x_10" 5 X-rays/ cm 2 . 


This spectrum is plotted in Figure 2, for various values of the electron 
e-folding energy E j. 
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MODEL II 

Finally, it should be noted that there are Jovian phenomena, such 
as the decametric radiation, which have no counterpart on earth. Expanding 
on the calculations of Warwick (6) , assuming that this is electron ,gyro- 

I 

radiation, suppose that the flux density observed at the earth due ;to 
decametric radiation is <j> keV/cm 2 sec Hz at a frequency f (Hz), and that 
the signal has a bandwidth Af (Hz) . If the radiation is beamed into a 

solid angle ft (sterad) , and the earth is a distance D (cm) away, the 

! 

total power radiated by the source is 


<p Af ft D 2 keV/sec . ' 

i 

A single electron of energy E (keV) radiates a total power 

4 x 10 -9 EB 2 sin 2 a keV/sec due to gyroradiation, where B is the magnetic 

field strength in gauss, and a is the angle between ”B and the electron 

velocity TT. In an electron spectrum electrons/cm 2 sec keV, the total 

at 

power radiated/unit volume is 

00 

4 x 10~ 9 B 2 sin 2 a J ~ dE keV/cm 3 . 

0 


dN 

T aking § - C 

(the nonrelativistic form for the velocity is justified again by the fact 
that the contribution from the exponential spectrum becomes very small for 
energies near me 2 , for the values of E q used here) , the total power radiated 
per unit volume becomes 


-E/F 

T e ' oJ electrons/cm 2 sec keV, and u = C 


cm/ sec 


4 x l(T 9 


B 2 sin 2 a 



E^ 3 / 2 T (3/2). 


c 
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The total volume into which the electrons precipitate is AjAR, where Aj 
is the area and AR the thickness of the precipitation region. Following 
Warwick (6), the thickness is found by assuming that the bandwidth of the 
decametric radiation is related to the change in the magnetic field as 
the electrons move closer to the planet 


AR = R 


AB 

3B 


Af 


Af 


3f " 3 x 2.8 x 10 b B 


where R is the planetary radius, AB is the change in the field strength, 
and the expression f * 2.8 x 10 6 B has been used for the gyrofrequency. 

t 

Thus 


or 


Aff2D 2 = 4 x 10" 9 B 2 sin 2 


C J 

a — \~2 


E 3/2 
oJ 


Af 


r t 3/2 > R 8.4 7l08g A . 


C ,E t 2 A t 
J oJ J 


r= <|>OD 2 

> E oJ RBsin 2 a 


4.5 x 1G 24 = 


1.4 x 10 32 



for Jupiter, using values in (6): 


(J> = 10" 19 watts/m 2 Hz = 6.24 x 10~ 8 keV/cm 2 

£2 = 8.63 x 10" 3 sterad (a cone of half-angle 3°) 

B = 8.9 gauss 

D = 4 A.U. = 6 x 10 13 cm 


sin 2 a has been taken to be 1, corresponding to a mirroring particle. 
From Appendix A," Formula TA, the X-ray spectrum at the earth is — ~ 


dn IL Chu) 
d (tiw) 


3.92 x 10' 3 



e -tlM/E oJ 

ha) 


X-rays/cm 2 sec keV 
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This is plotted in Figure 3, for various values of the e-folding energy 
E .• 

°J 

X-RAYS FROM Io 

Goldreich and Lynden-Bell (7) have estimated that, since Io is immersed 
in the Jovian radiation belts, it should be a source of kilovolt X-rays; 
they give a flux of roughly 10" 7 X-rays/cm 2 sec at the earth. Thus this 
is comparable to the flux of Model IA. ! 

SATURNIAN X-RAYS 

i 

There have been reports (8, 9, 10, 11) of sporadic non-thermal 
radio emission from Saturn; theoretical estimates of Saturn's magnetic 
moment based on these observations (12, 13, 14, 15), or on general 
scaling relationships (16, 17) vary widely. For the following section, 
an intermediate value of M g = 2 x 10 29 gauss-cm 3 (or 2.5 x 10 3 M £ ) is used. 

This gives a surface field of about 1 gauss at the equator, and a magneto- 
spheric radius of 39 Rg at the subsolar point. As the radius of Saturn 
is r = 9.45 r p , the area of the precipitating region should be 2.68 x 10 17 
cm 2 . 

MODEL IA (SATURN) 

Following the calculations for the Jovian Model IA, the scaling 
factor relating Saturnian to terrestrial intercepted energy becomes 

2/3 

= 9.2, so AgCgE^g 2 = 9.2 A^CpE^ 2 = 1.94 x 10 26 keV/sec 


m s r e 

m e r s 2 
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As the distance between Saturn and the earth is 10.5 astronomical units 
(1.57 x 10 14 cm) at MJS launch, the Saturnian X-ray spectrum is about 
one order of magnitude below the Jovian spectrum, or, from Formula IA, 
Appendix A, 


dn(fiw) _4u/E 0 s 

— = 1.5 x 10" 9 - 


d(fiio) 


-ho) 


X-rays /cm 2 sec keV 


This is plotted in Figure 4, for various values of the e-folding energy 

i 

F ! 

oS* 

MODEL IB (SATURN) j 

i 

Following the calculation for the Jovian Model IB, we obtain a scaling 

i 

factor 

I M s r E 3 V 

— 3 = 8.77 

\ r s 3 M E 


so 


C c E _ 3 / 2 = 8.77 CpE c 3 / 2 
S oS E ofc 


or 


C S E oS 2A S = 8 - 7? C E E oE 3/2 ^oS - 5 - 2 ° X 1027 i 

! 

From Appendix A, Formula IA, the X-ray spectrum at the earth is 1 


dn IB (fioj) 

dfruj) 


= 4.03 x 10- 


“J^oS 


; -ho)/E 0 s 
" 1tu3 


X-rays/cm 2 sec keV 


at launch. This is plotted in Figure 5, for various values of the e -folding 
energy E qS . 




DNIB/DHW (PHDTQh~/'Cf12 SEC KEU) 
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KLtE0*lO>« MWtOHOlf 0*90) » *8<JWtS(C0*SO)» M.»«<CO*iOO) 
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MODEL II (SATURN) 

i 

I 

This is based on the same argument as the Jovian Model II. From 
(18), the value = 10" 22 watts/m 2 Hz = 6.24 x 10" 11 keV/cm 2 is used. 

The value of B is taken to be 2 gauss near the surface of the planet 
at the poles, and, as a guess, the same value of is used for Jupiter. 
Thus 


c s E os 2A s - 4 - 5 * 1 ° 24 N p7s 


= 5 X 10 


30 


oS ' 


From Appendix A, Formula IA, the X-ray spectrum at the earth is 


dn (ho)) — -ha)/E oS 

-SW-' 3 - 87 * 10-5 J*os Hsr- 


X-rays/cm 2 sec 


This is plotted in Figure 6, for various values of E Q g. 


keV. 



I 
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II. INSTRUMENTATION FOR A FLYBY X-RAY EXPERIMENT 

Much of the work done for this phase was written up in "Instrument 
Report: Planetary X-ray Experiment"; we have therefore included this 

report as Appendix B. As it seemed desirable to measure not only brems- 
strahlung X-rays, but also fluorescent X-rays from the satellites of 
the outer planets (and therefore obtain composition information) , we 
considered the possibility of using very high resolution Si (Li) detectors. 
Two serious problems were found, however: 

1 a) In the Grand Tour configuration, radiation damage from both 
RTG's and planetary trapped radiation would have decreased 
the useful lifetime of the detectors to less than the mission 
lifetime. j 

b) The background counting rate from the RTG's would be very high. 

We experimented with Si (Li) detectors and fast rise time pre- 
amplifiers, attempting to do pulse shape discrimination, which 
would help to solve the background problem. We concluded that 
this technique was not feasible, and turned our attention to 
proportional counters. On the basis of the upper limit calcu- 
lations in the Instrument Report, we determined that a propor- 
tional counter was capable of operating reliably, and with a 
reasonable background rate, on a Grand Tour type spacecraft. 

— -One point which required further investigation ^was the problem of 
Compton electron production in the proportional counter walls by RTG gamma 
rays. We developed a Monte Carlo program to study this effect, and to 
improve on the upper limit calculations in the Instrument Report. 
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Figure 7 
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At the same time, we began construction of a prototype detector; 
this is shown in Figure 7. It consists of a proportional counter with 
a passive collimator restricting the field of view to about 7° x 7° FWHM 
and a modulation collimator consisting of 2 electrodeposited tungsten 
meshes. The purpose of the modulation collimator is to provide X-ray 
source size information down to about 0.1°. We tested this double 
collimator with various size light sources; the results of the test are 
shown in Figure 8. As the collimator is rotated, the detector response 
is a roughly triangular envelope (the response of the 7° x 7° collimator) 
containing a series of finer triangles (the modulation collimator response) . 
The peak to valley ratio of the modulation collimator response indicates 

i 

the source size. 

Finally, we tested the prototype detector at JPL, in the presence of 
the simulated RTG; the results of this test agreed well with the background 
estimates we had made using the RTG description in the MJS AFO package, in 
conjunction with our Monte Carlo results. 

The detector which we propose to fly on the MJS mission is shown in 
Figure 9; it is divided into two separate gas chambers for redundancy, 
and each gas chamber has smaller anticoincidence chambers around the 
sides to eliminate the Compton electron contribution. We have estimated 
the total background count rate in this detector to be 2.2 counts/sec on 
the MJS spacecraft. 



INTENSITY INTENSITY INTENSITY 





Figure 9 



III. STATISTICS AND METHODS OF OBSERVING JUPITER X-RAYS 


Let n be the number of counts/second received by a given detector 
in some energy range due to an X-ray source (formula 2A of Appendix A) , 
and let B be the number of counts/sec, due to background in the same 
energy range. In a time T, it is desired to measure both n and B, and 
to derive n with the maximum accuracy. The count rate is thus n + B 
during the source observation, and B during the background observation. 
Suppose a fraction a of the time T is devoted to the source, so that 


(n + B) a T ± [(n + B) a T] 1 / 2 

counts are received during this observation period. Then a time (1 - a 
is left for the background observation, and 

B (1 - a)T ± [B Cl~ oOT] 1 / 2 


counts are received during that period. The number of source counts is 
derived by multiplying the number of background counts by ^ and 

subtracting, to get 


n a T ± ! (n + B) a T + 


BT 


1/2 


The maximum accuracy of this quantity is obtained for the minimum 
value of the quantity 


f n + B > a T + 


~cr 


BT 


1/2 


jn ct T 


and this is obtained for 
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Requiring the number of source counts to be 3 standard deviations 
above zero (for 99% confidence) means that 



n + B 1_ 1 ) 

a + n 2 1 - a | 


We will consider only the case n << B, to derive an upper limit 
to the time needed. For this case a « 1 - a « 1/2, i.e., equal times 
are needed for source and background observations; the total time needed 
is 

T = 36 - 2 

n 2 

Alternatively, in a given time T, the minimum detectable count rate is 



If the detector area is A d , and the average efficiency over the energy 
range of interest is e, the minimum detectable flux is 



Using these formulas, we have considered many possible methods for 
observing the models of Jovian X-ray emission discussed in Section I. The 
most sensitive of these are the experiments which operate outside the earth's 
atmosphere, in the energy range below about 30 keV; three of them are 
explained below, and summarized in Table I. 

SATE L L~ITE~PR0P0RT I ONA L“C0UNTE r— 


This is the AS§E experiment on the SAS-A (UHURU) satellite (19, 20). 
The satellite is spin stabilized, with the spin axis reoriented every 
day; the instrument is collimated to 5° x 5° FWHM, so the duty cycle 
for observing any given source is 10°/360° = 0.028. The response pattern 



Table 
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of the collimator is triangular, so the total observing time must be multi- 
plied by 1/2 to obtain the effective time. As the earth occults any 
given source for a total of 12 hrs/day as seen from the satellite, the 
effective time available on a source is (12 hrs) (0. 028) (1/2) = 0.168 hrs, 
or 10 minutes. The area-efficiency product is £840 cm 2 x 10% = 84 cm 2 . 

The actual efficiency function e(fi<i)) used in Formula 2A is given in the 
literature. 

SATELLITE GRAZING INCIDENCE TELESCOPE 

This experiment is the AS$E grazing incidence telescope to be 
flown on HEAO-C (21). It is capable of 2 arc second resolution over 
a 1° field of view; since the disc of Jupiter subtends about 50 arc 
seconds, the telescope could operate in two different modes. First, a 
mapping mode, in which the disc is scanned one resolution element at a 
time, until all (50/2) 2 = 625 elements have been scanned; second, simply 
an X-ray detection mode, where all 625 elements are scanned at once. 

In the mapping mode [referred to as the point source mode in (21)], 
each resolution element gives 50 background counts in an observing period 
of 5 x 10 4 seconds (14 hours), for a background counting rate of 10“ 3 
couhts/second/resolution element. Using equal observing times for Jupiter 
and background, the total observing time must be divided into 2 x 625 = 1250 
equal segments, or 40 seconds/segment, during which a resolution element 
will contribute 0.04 counts. If we require the total probability of a 

random-error-simuTating-a-Jupiter - X-ray^event-to be O.OT for-the-625 

Jupiter observations, Poisson statistics may be used to find the required 
number of counts in any Jupiter segment; for a mean of 0.04, three counts 
are needed for 1.6 x 10" 5 probability (= 0.01/625) of random simulation. 
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The sensitivity derived here for the mapping mode is much less than 
the sensitivity given in (21), where it was assumed that the background 
was well known, and that a total of 5 x 10 4 secs would be devoted to a 
single resolution element. This type of observation would be used if, 
for example, it were desired to correlate an X-ray source with an 
optical source whose position were well known. 

In the detection mode, the total observing time would be divided 
into just two segments - Jupiter and background. The background 
counting rate would be 0.625 counts/second. Dividing the time equally 
between background and Jupiter, the total number of counts in a segment 
would be 1.5 x 10 4 , so Gaussian statistics may be used, and the formulas 
derived at the beginning of this section may be used. 

MJS SPACECRAFT PROPORTIONAL COUNTER 

This is a Xenon-filled counter, as described in Section III and 
Appendix B, with an area-efficiency product of 15 cm 2 . It has a modulation 
collimator which makes possible the mapping of a planet in X-rays to about 
1/10°; at entry into the Jovian magnetosphere, then, its resolution is 
comparable to the grazing incidence telescope operating in earth orbit. 

Two cases are considered: the sensitivity of the counter at launch, 

5.5 A.U. from Jupiter, and the sensitivity just before entry into the 
magnetosphere, 0.02 A.U. from Jupiter. So that this experiment can be 
compared to the grazing incidence telescope, the observing time has been 
taken to be 14 hours, for the 5.5 A.U. case. At 0.02 A.U., all models 
become detectable, and the only question is the amount of time needed; 
this is tabulated in the last three columns of Table I. 
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CONCLUSION 

As the table shows, the most sensitive earth-based detector is the 
HEAO grazing incidence telescope; by comparison, the MJS proportional 
counter is not particularly sensitive if used near the earth. The over- 
whelming advantage to the MJS experiment is its ability to get very 

close to Jupiter, where the X-ray fluxes can be observed in short 

» 

integration times. If simply the detection of Jovian X-rays were 
desired, this would be theoretically possible for a number of near earth 
experiments; but if time variations, detailed mapping, and both dayside 
and nightside observations are required, it is absolutely essential to 
use a flyby experiment. 
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APPENDIX A: PRODUCTION OF X-RAYS BY BREMSSTRAHLUNG 

A single electron of energy E, traversing a thick target of atomic 
number Z, produces a bremsstrahlung X-ray spectrum given by 


T©- >- 2 x10 ' 6 2 (fit - 


photons/keV 


where dn(fiu)) is the , number of photons of energy h to "ttio + d(fia). This 
result, known as the Kramers formula, is derived in detail in (1); it is 
valid in the nonrelativistic energy region, and it has been applied suc- 
cessfully to the problem of X-ray production in the earth's auroral 

zones (22). If the number of electrons with energies between E and 
dN 

E + dE is -tft dE, the number of photons/keV at an energy -hm becomes 
at 


dn(fico) 

d0to) 


1.2 x IQ" 6 Z 


/ 


E 

dial" 


dN 

dE 


dE 


An energy spectrum for precipitating electrons which is frequently 
observed in the earth's auroral zones is 

^ = C e” E ^ E ° electrons/cm 2 sec keV . 
dE 

The parameters C and E q (the e- folding energy) are often observed to 
vary on a time scale of seconds or minutes; however, a typical spectrum 
can be taken to be 


4rr = 7 x. io 7 e electrons/cm 2 sec keV 

dE — : 
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Although other functional forms can be found which fit the experimental 
data (22), the exponential spectrum will be used for the following 
calculations. This gives a photon spectrum at the production region 


dr ) = 1.20 x 10 -6 ZCE 2 
d (iSw) o 


3 --Ko)/E 0 


X-rays/cm 2 sec keV 


To find the X-ray spectrum a distance D away, this expression must be 
multiplied by A/4?rD 2 , where A is the area into which the electrons pre- 
cipitate. The assumption implicit in this factor is that the X-ray 
emission is isotropic; this is justified, not by the radiation pattern 
of the electron, which becomes strongly peaked in the forward direction 
at higher energies, but rather, by the fact that electron scattering 
renders the electron flux isotropic in a distance which is small compared 
to the electron range (1). Thus the X-ray spectrum becomes 

. 2.4 X 10-« CE o 2 X-rays/cm 2 sec keV (1A) 

adopting the value Z = 2. Values of C, E q , and A are derived in Section I. 

Another quantity which is of interest is the total number of counts/ 
second which this spectrum will yield in a detector. This is obtained by 
multiplying the spectrum by the detector area A^, and the detector effi- 
ciency e(ftw). Thus the total number of counts/sec in a detector, in the 
energy window "Rte i -* > is 

d(Ru>) (2A) 


nfifcoi -* ft 0)2) 
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INSTRUMENT REPORT 
PLANETARY X-RAY EXPERIMENT 
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PART I 


INTRODUCTION 



The basic scientific motivation for an x-ray experiment is to 
discover the dynamics of the Jovian and Saturnian magnetospheres. 

Since the Earth's magnetosphere is now reasonably well known, 
important comparisons can be made for three planetary objects. 

The x-ray experiment is intended to study particles and plasmas in 
the magnetospheres and the processes that act on them (e. g. , waves, 
convection). During the encounter with the magnetosphere a profile, 
essentially at one instant of time, is made by the direct particle 
counters. There are three major difficulties with such measurements: 

1. The spacecraft moves so rapidly that time variations of 
plasmas, particles and fields cannot be determined except 
over a narrow range of frequencies. 

2. The trajectory cannot explore all regions of interest. 

3. It is probably not possible to measure precipitating particle 
fluxes. These particles carry information about magneto- 
spheric phenomena (waves, auroras, substorms). 

The x-ray experiment is intended to remedy these difficulties. The 
hope is that x-ray fluxes could be detected even months before 
encounter, then measured from that time until encounter and then for 
months after the encounter. In that way the pattern of magnetospheric 
activity could be determined. 

The x-ray experiment is one of a group of experiments that can 
contribute knowledge of planetary magnetospheres. The other experi- 
ments are those to measure fields, particles and plasmas while in 
the magnetosphere; and radio noise emission measured inside as well 
as outside the magnetosphere. A secondary scientific goal is to make 
use of the long interplanetary cruise times to study long term variations 
and transient events in a few cosmic x-ray sources. 
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Design of a Planetary X-ray Experiment 

Some of the design constraints that an outer planetary x-ray 
experiment must satisfy: 

1. Reliable operation for periods up to 10 years. 

This requirement demands that gain drifts be very low 
and that they be corrected in-flight. 

2. For very long missions using MHW RTG power sources, 
the detectors must be insensitive to radiation damage 

by neutrons and protons. Furthermore, the neutrons and 
gamma rays from the RTG sources must not contribute 
very much to the background rate of the detector. 

3. The weight should be limited to 3 kg, although more 
sensitivity may be gained from larger area detectors. 

4. The detectors should have good energy (~30%) resolution 
down to ~2 keV photon energy. 

5. The angular resolution should be about 0. 1° in order to 
examine satellites of the major planets, and to find the 
pattern of the x-ray emission on the planet. 

The following types of instruments were considered: 

1. Grazing incidence focusing optics 

a. Geiger-Mueller tube detector 

b. Channeltron detector (for x-ray energies below 1 keV) 

c. Solid-state detector (small area) 

d. Proportional counter (small area) 

2. Scintillation counter with passive collimator and modulation 
collimator 
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3. Solid state detector array with passive and modulation 
collimators 

4. Proportional counter with passive and modulation 
collimators 

a. Single wire type 

b. Multiwire for background rejection 

Channeltrons were ruled out because of their short lifetimes 

(10^ counts) and large gain drifts which are difficult to correct for 
beyond a certain point. The scintillation counter, although feasible 
from most points of view, does not have the energy resolution of 
the proportional counter. Also, the weight per square centimeter of 
sensitive area is more favorable for other detectors. The Si(Li) 
detector appeared promising at first: its excellent energy resolution 

might allow surface chemistry experiments on the outer ‘planet 
satellites without atmospheres. 

The grazing incidence x-ray telescope (for example, the 
hyperboloid-paraboloid type) would allow relatively large x-ray collec- 
ting areas to be achieved with small area detectors. In addition, 
the small field of view would make it possible to do rather detailed 
mapping of x-ray emission regions. However, at the present state of 
the art, this type of collimator does not have a high efficiency for 
x-ray reflection; conventional collimators give equal effective areas 
for less weight. Therefore, a wire grid modulation collimator, with 
resolution of about 1/10°, appears to be a better choice. In addition, 
rough collimation, either active or passive, down to a field of about 
5° square, would be required. 
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With the choice thus limited to two detectors, and wire grid 
+ active collimator or wire grid + passive collimator, two design 
studies were undertaken: one for a solid state detector, and one for 

a proportional counter. These are included as Parts II and III. 

It was recognized from the start that the MHW RTGs would be a major 
problem for the x-ray experiment; therefore, the design studies con- 
centrated on the question of background counting rates, and how to 
reduce them to an acceptable level with a minimum of weight. At 
the same time, because the spacecraft design was not finalized, and 
the MHW RTG radiation levels not known accurately, the emphasis 
was placed on deriving upper limits to MHW RTG -induced counting 
rates, and identifying problem areas, rather than obtaining accurate 
estimates of counting rates. This must await the final design of the 
spacecraft and RTGs. 

The design study for the solid state detector showed that the 
background counting rate would be quite high, and that a significant 
portion of it would be caused by RTG neutrons. In addition, radia- 
tion damage due to RTG neutrons, as well as radiation belt protons, 
would be severe with currently produced detectors. Since radiation 
hardening techniques are being developed, there is at least some 
hope that the radiation damage problem can be overcome. Therefore, 
it was decided to proceed experimentally with solid state detectors, 
to determine whether rise time discrimination could be used to solve 
the background problem associated with the RTG neutrons. We 
have tested a Si(Li) detector with two fast rise time (50 ns) preampli- 
fiers; one of these was constructed in our laboratory, and the other 
is on loan from General Electric, Daytona Beach. Initial results are 
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not promising: it appears to be impossible to do rise time discri- 

mination in the energy region where the x-ray experiment must work 
(below 30 keV), due to the very short length of the ionization tracks. 

The proportional counter design study is more encouraging; the 
major problem area was identified as due to RTG gamma rays. 

Direct interactions in the counter gas, as well as interactions in the 
counter body giving rise to electrons which deposit energy in the 
counter are both important effects. However, the use of multiwire 
counters and/or rise time discrimination can solve these problems; 
both techniques are currently in use in the energy range of interest. 
Some passive shielding from the RTGs will be provided by the space- 
craft itself; further shielding, either active or passive, will be 
required at the detector. Since the spacecraft design is still some- 
what uncertain, an accurate evaluation of the gamma flux at the 
detector cannot be made, and the amount of additional shielding 
needed is not known. On the basis of upper limit estimates, however, 
the design study shows that a shield weighing about one pound will 
reduce the background to an acceptable level. This shielding is not 
required for RTGs as they would be operated on MJS missions. In 
order to make a more accurate evaluation of active versus passive 
shields, comparison of geometries, etc. we are currently developing 
a Monte Carlo program. 

Experimentally, two Xenon filled proportional counters are 
being tested; a third, Krypton counter will be delivered Shortly. ~ 

A delay-line shaping amplifier is being used in conjunction with a 
pulse shape analyzer to do rise time discrimination. A passive, 
slat collimator with a built in tungsten wire modulation collimator is 



now being designed; this will be constructed in our machine shop 
and tested in the near future. 



PART II 


DESIGN STUDY FOR AN X-RAY EXPERIMENT 
USING SOLID STATE DETECTORS 
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I. THE DETECTOR 

Figure 1 shows the general configuration of the detector assembly- 

used in this study. Four Si(Li) detectors, each .2 cm thick and having 
2 

.5 cm area, are used for the 1.5-10 keV range; two Si(Li) detectors, 

2 

each .2 cm thick and having 5 cm area, are used for the 10-30 keV 
range. 

These six detectors are placed in the bottom of a square, well -type 
collimator measuring 10 X 10 X 15 cm, to give approximately a 37° 
(FWHM) field of view. In the well collimator, above the detectors, is 
a tungsten wire modulation collimator, with wire spacing approximately 
equal to wire diameter. (Performance of the modulation collimator, 
selection of wire diameter and spacing, etc., will be discussed in a 
later section. ) This rather simple geometry has been chosen to simplify 
the calculations of the background counting rates. Other detector con- 
figurations should give approximately the same results, as long as the 
surface areas and fields of view are approximately the same. 

The power source for the Grand Tour spacecraft, a 2200 W(Th) 

RTG, has been assumed to be 10 feet away from the x-ray experiment, 

with nothing in between the two. The neutron and gamma radiation levels 

2 

are thus obtained by 1/r extrapolation of the levels given by Trainor 1 
and Noen and Anno. 2 

All interactions in the collimator (if it is active) are assumed to 
give rise to an anticoincidence signal. Since such a signal will contri- 
bute to the detector dead time, the background has two effects, First, 
interactions in the detectors, if they are in the 1.5-30 keV range, 
contribute to the detector background counting rate; second, shield 
interactions contribute to the detector dead time. The relative importance 
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of these two effects is as follows. Suppose that an x-ray source, whose 
true strength is S counts /sec, must be detected in the presence of a 
background whose true strength is B counts/sec. The detector dead 
time is bt, where b is the number of interactions /second producing an 
anticoincidence pulse, and t is the length of the gate signal for such a 
pulse. Then in a time interval T, ST-STbt counts are collected from 
the source, and BT-BTbt from the background. The source will be 
detected if ST(1 - bt) £ 3a bkgd = 3/BT(l - bt) . Thus T 2 9B/(S 2 (1 - bt)], 
and the factor B/(l - bt) determines the time needed to detect a source. 

II. BACKGROUND COUNTING RATES THROUGH FORWARD APERTURE 
A. Diffuse Cosmic X-ray Flux 
The spectrum of the diffuse flux is: 3 

~ = 10 E ~ 1 * 5 / cm 2 sec sr keV, 1 < E < 10 keV 

d.ii« 

= . 16. 8 .E" 1, 75 /'cm 2 sec sr keV, 10 < E < 40 keV 
= 128 E -2 ’ 3 /cm 2 sec sr keV, 40 < E < 1000 keV 

The total fluxes in energy ranges of interest to this experiment are: 

10 2 

N(l. 5-10) = J jg dE = 10/cm sec sr 
1.5 

: ~ ~ 7 ' "• “T - ■■ 

N(10-30) = 2.5/cm sec sr 

N(>30) = 1.18/cm 2 sec sr 

2 

= . 098/cm sec sr 


N( >200) 



10 . , 

2 

Using geometrical factors of . 052 cm sr for 1.5-10 keV, and .26 
2 

cm sr for 10-30 keV, which assume 100% efficiency for the shield 

and modulation collimator, the counting rate in the 1.5-10 keV 

detector is . 52 counts /second. In the 1.5-10 keV range, the probability 

2 

of an interaction in . 2 cm (.47 g/cm ) Si is close to 100%. Above 

15 keY, however, the efficiency begins to fall rapidly, and the spectrum 

must be weighted by 1 -exp(-p^,x), where pvp = total attenuation 

2 

coefficient, and x = .47 g/cm . 

30 , / -ij, t X\ _ 

J dE - e J dE ~ 2/cm sec sr 

10 

Thus the count rate is 2(. 26) = . 52 counts /second. 

The diffuse flux can also contribute to the background count rate 
by scattering and leaving 1.5-10, or 10-30 keV in the detectors. The 
graph of Figure 2 shows for . 2 cm thick Si: 

1, The probability of a photon interaction of any kind, 1-e 

2. The probability of a Compton interaction leaving 1.5-10 keV, 


a Compton, 1,5-10 
a total 



where the cross sections 


1 0 keV 


Compto 


n, 1.5-10 I 


d a 


1.5 keV 


dE 


dE 


a , = total cross section for interaction of any kind 
total 
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3. The probability of a Compton interaction leaving 10-30 keV, 


Compton, 10-30 
a total 


0 - 


and Q^, are given by Grodstein 4 ; M , (^ om p^- on is given by Nelms. 5 

Figure 3 shows the diffuse x-ray spectrum, and the product of 

this spectrum with the three functions of Figure 2. Integrating the 

bottom two curves of Figure 3, the contribution to the 1.5-10 keV count 
- 2 2 

rate is 7 x 10 /cm sec sr, and the contribution to the 10-30 keV 

-2 2 

count rate is 2.5 X 10 /cm sec sr. Using geometrical factors of 
2 2 

.208 cm sr and 1.04 cm sr for 1.5-10 and 10-30 respectively 
(derived by assuming that the modulation collimator is transparent, 
while the shield is still 100% efficient) the count rates are 
.0145 counts/sec, 1.5-10 keV 
. 0256 counts/sec, 10-30 keV 


While these contributions could be reduced by the use of active 
anticoincidence (since the scattered photons would have a chance of 
interacting in the bottom of the well and being detected) they are already 
small enough to be neglected. 


Table 1 

Background Contributions Due to Diffuse X-rays 
(Passive Collimator) 

• ' ' • *' T. 5-10 keV ~ 10-30 keV 

Photoelectric interactions .52 counts/sec .52 counts/sec 

Compton scatter .0145 counts/sec . 026 counts/sec 
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III. LEAKAGE FLUXES 

This section deals with fluxes originating outside the field of view 
of the detectors; i. e„ , fluxes which leak through the shield. 

A. RTG Gamma Rays 

Gamma and neutron fluxes from the RTG are given in Trainor 1 

and Noen and Anno 2 ; the gamma spectrum 2 for 0° and 3 feet is shown 

2 

in Figure 4 extrapolated by 1/r to 1 0 feet. 

The background counting rate in the detectors due to this flux 

-)J*jiX / cos 9 . 

will be proportional to A cos 9 (1 - e 



where A is the detector area, x the thickness, the total absorption 

coefficient, and 9 is the angle between the incident flux and the detector 

-g^,x/cos 9 

normal. As long as (p r j,x/cos8) ^ . 2, 1 - e ~ gx/cos9 , to 

15% or better. Thus the counting rate becomes proportional to Agv-pX, 

and does not vary with the orientation of the detector. The range of 

■validity of this approximation, as a function of incident flux energy, is 

_ x ' -p T x/ cos 8 

shown in Figure 5. Also, since e s 1 -x, Ag^,x s A cos 9(1 - e ), 

-so-t-he-use- of the— approximation— will- give— an upper limit— to^ the counting 

rafe even when it is not strictly valid. 
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Thus Figure 4 also shows the product of the RTG gamma spectrum 
by the probabilities of Figure 2 and the areas of the 1.5-10 and 10-30 
keV detectors. These two graphs show the spectra which would be 
present in the detectors if there were no collimation at all; the total 
counting rates are 8 counts /second, 1.5-10, and 49 counts /second, 10-30. 
The spectra of Figure 4 show that the main contributions are from photons 
with energies of 10-200 keV. In this energy range, either an active or 
a passive collimator could be used. An active collimator can attenuate 
photons either by photoelectric or Compton interactions, and the photon 


can interact either before or after depositing a count in the Si detector. 



Thus the intensity I - I Q e ^ X » an< ^ the attenuation factor ^ ac ^ ve 

— 

A passive collimator has only one chance to stop an incoming 
photon; before it reaches the main detectors, and the interaction must 
be photoelectric. Thus the attenuation factor 
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£ 

passive 





where Up^. is the attenuation factor for photoelectric interactions only. 

The weight W of a thin wall, well-type shield (Figure 1) is 
W = p V, where p is the density and V is the volume. If A is the surface 
area and t is the thickness in cm, W = p At. But p t = x, the thickness 
in g/cm , so W = Ax. 

For an active shield, 

x = - ^ "A - , so W = - ~ ln(l-f). 

2 [x 2 (1 ' ' 


For a passive shield. 

In i 1 


M*PE 


so W = - 




In ^1 


^VE 


Examples 

1. A Silicon shield which actively stops 80% of 200 keV photons: 

W — — ln(l-. 80) = 6 . 35 g/cm 2 A 

2 x . 127 cm g 

2 

Thus the shield must be 6.35 g/cm thick; using an area 


A = 700 cm 2 , W = 4400 gm 

=10 pounds. 


2. A Csl shield which actively 

stops 80% of 200 keV photons: 


W 7 — 

ln(l-.80) = 2..9A 


2 X . 273 cm /g 
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This shield would therefore be 2. 9 g/cm thick, and would weigh 
about 2100 gm, or 4.5 pounds. 

A lead shield which passively stops 80% of 200 keV photons: 


W = 


A 

. 942 crn^g 


In ( 1 


.942 
. 80 



This shield would be about 3 g/cm thick and weigh about 

2100 gm, or 4. 5 pounds. 

The above figures represent upper limits, since the incoming 
photons were assumed to be normal to the shield. It would appear that 
the lead shield is preferable, since it does not require any supporting 
electronics which would add to the weight and complexity of the system. 
However, the fluorescent x-rays produced in lead and Csl must be 
compared, since they will contribute to the background counting rate. 

An upper limit is now derived for this contribution. Assume: 

1. That every photoelectric interaction in the shield gives 
rise to a K x-ray (i. e. , ignore the less energetic L, M, 
etc. x-rays, and the Auger effect). 

2. That every fluorescent x-ray leaves the atom in the 
direction which gives it the shortest escape path through 
the shield. 


x ► 
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The following quantity is now computed: 

{Probability of an incoming x-ray of energy E going a distance 
y with no photoelectric interaction) X. 

(Probability of x-ray of energy E photoelectrically interacting 
in a layer dy thick) x 

(Probability of a fluorescent x-ray of energy E' going a distance 
x-y with no photoelectric interaction). 

For any x-ray, Probability of Interaction + Probability of No 
Interaction = 1. And Probability of Interaction = Probability of Photo- 
electron Interaction + Probability of Compton Interaction + Probability 
of Rayleigh Interaction + (higher order cross terms). Using P(I), 

P(NI), P(PE), P(C), and P(R) for the above probabilities, and neglecting 
the cross terms: 

P(PE) + P(C) + P(R) + P(NI) = 1, or 

P(PE) + P(no PE) = 1, where 

P(no PE) = probability of no photoelectron interaction 
= P(C) + P(R) + P(NI) . 


Since 

P(NI) 

= e-»*y 

P(I) = 1 - e ^ y , and 


P(C) 

II 

- e-** y ) 


P(R) 

= — (1 

p 

- e^ y ) . 


P(PE) 

^PE 
. (1 

P 

- e^ y ) 


where p = + p^, + Mp is the total interaction coefficient, and 

, p^., p^ are the coefficients for photoelectric, Compton, and 
Rayleigh interactions. 




Thus P(no PE) may be written as 


P(no PE) = ■ — E (1 - e"^ 7 ) + e~^ y 






The probability of a photoelectric interaction in a thickness dy is 


^PE -pdy x j 

— ( - C > = I'PE dy 


The probability of no photoelectric interaction in a distance x-y, for 
an energy E' (coefficients p'^) is 


^PE “P'(x-y) 


PE 


M- 


Multiplying the above three probabilities together, and integrating over 
y values from x to 0, we get for the probability that an incoming photon 
of energy E produces a fluorescent photon of energy E' which escapes, 


f ^PE 

hPE 1 p 


n* 


PE 


-11X -p'x 
e - e 

M.' - M- 


VJt- 


+ 


PE 


W - (J- 


PE 


1 - 


e"“ x \ + 
/ 


» - ^PE ^ : ^PE - ^PE ^-PE 1 - e Vx 1 

p jj, 1 X p ~p7 p 1 J 

For lead, take E = 200 keV, E 1 - 80 keV, x = 3 g7cm^, ppp = .80 crri^ 

2 2 2 
p = .94 cm /g, p'pg = 1.29 cm /g, p* = 1.66 cm /g. The probability 


then becomes . 32. 
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For Csl, take E = 200 keV, E' = 30 keV; x = 2.9 g/cm 2 , 

Mpp =.250 cm 2 /g, p = .272, M-'p^ = 6.33, p, 1 = 6.34. The probability 
for Csl becomes .02. 

How many of these escaping photons interact in the detectors? 

To get an upper limit to this number, assume that the escaping photons 
are distributed isotropically, and let n be the number which escape 
per cm 2 per second. 



51PE VIEW 

A thin strip sdl emits nsdl photons /second, and the flux at a distance 

2 2 22 
(1 + r r from this strip is nsdl/4n(l +r ). The detector (area A) 

presents an effective area A coscp to this flux, and the probability that 

-UX/ COS CD 

a photon will interact in the detector is 1 - e . Thus as an 

xisd.1 

upper limit to the contribution from the thin strip, we have Apx. 

Integrating over 1 values from 0 to L, this becomes 


nsAux 1_ 
4tt r 


tan 


L 

r 


Typical values are r = 5 cm, . L = 15 cm and s = 10 cm, so we get 
.2n A(i x, 
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From Figure 4, the incident flux of 200 keV photons is 

2 2 7 

1.5 X 10 /cm sec. Thus for a lead shield, n = (. 32)(1. 5 X 1 0 ) = 

2 

48/cm sec. From Figure 2, the contribution to the low energy 

_2 

detectors is (48) (2) (. 2)(3. 5 X 1 0 ) = .67 counts/sec and the contribution 

-2 

to the high energy detectors is (48) (1 0)(. 2) (3 X 10 ) = 2.9 counts/sec. 

2 2 

For the Csl shield, n = (. 02)(1.5 X 10 ) = 3 /cm sec and from 
Figure 2 the low energy contribution is (3) (2) (. 2)(3 x 10 - ^) = .036 
counts/sec. Since a 30 keV photon can leave at most several keV to an 
electron in a Compton interaction, the contribution to the high energy 
detectors will be through photoelectric interactions only. From Figure 2, 
the probability of an interaction of any kind is 50% at 30 keV; the ratio 
of photoelectric to total cross-section is 1.17/1. 44. 4 Thus the contribu- 
tion to the high energy detectors is (3) ( 1 0)(. 2)(. 5) (1 . 17/1 . 44) = 2.4 
counts/sec. However, most of the Csl fluorescent x-rays will be gated 
out by means of the initial interaction in the shield. 

Since 80% of the 200 keV photons interact in the shield, the contri 
button to the dead time is (. 8) ( 1 . 5 X 1 0^) (1.5 X 10^ cm^) = 1.8 X 10^ 
counts/sec. 

The two shields considered above both stop 80% of the incident 

200 keV photons; to complete the comparison, the leakage at other 

energies must be computed. Thus the graph of Figure 6 shows the 

2 

probability of no photoelectric interaction in 3 g/cm of lead, as derived 

on page 14, and the graph of Figure 7 shows the probability of no 

2 

interaction at all in 2 X 2. 9 g/cm of Csl, as defined on page 21. 

Finally, the product of these probabilities with the 1.5-10 keV graph of 
Figure 4 gives the flux of photons which leaks through the shield and 
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scatters in the low energy detector, leaving 1.5-10 keV. This is shown 
in Figure 8, where a direct comparison of the active and passive shield 
counting rates (excluding the effect of fluorescent x-rays) may be made. 
Since the contribution from 100-200 keV photons is about an order of 
magnitude larger than any other contribution, and the shields have 
identical stopping powers in the 100-200 keV range, the difference 
between the two is negligible. The comparison is summarized in Table II. 

Table II 

RTG Photon Count Rates 

Passive Shield Active Shield 

(keV) (keV) 

1. 5-10 10-30 1. 5-10 10-30 

Leakage 1.6 c/s 9 c/s 1.6 c/s 9 c/s 

Fluorescence .67 c/s 2.9 c/s 0 0 

B. RTG Neutrons: Direct Effects in Si Detectors 

2 

The RTG neutron spectrum of Figure 9 was obtained by 1/r 

extrapolation to 10 feet of the spectrum given in Trainor. 1 It will 

be shown in section C that both the active and passive shields are 

almost transparent to neutrons; thus the entire spectrum of Figure 9 

may be considered to be incident on the Si detectors. The possible 

interactions of neutrons with Si are: 

28 28 

1. Si (n, n) Si (elastic scatter) 

28 28* 

2. Si (n,n') Si (inelastic scatter) 

3. Si^ (n, p) Al^; threshold energy 3.86 MeV. 6 

Al - 8 T Si , T^ =2.3 minutes. 


Dead Time 
1. 8 X 10 4 c/s 
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4. Si^ (n,#) Mg“^; threshold energy 2.66 MeV 6 

Figure 10 shows the cross-sections for these interactions as a 

function of energy 7 and Figure 11 shows the probabilities of these 

2 

interactions in .47 g/cm Si (E. g. , the probability of an elastic scatter 

is (a . ■ ' ,, /a) (1 - e ax ) where x is the thickness of the material, 

elastic scatter ■ 

a is the total cross-section for an interaction of any kind, and 

O', ■ .. is the eleastic scattering cross-section.) 

elastic scatter 



3 

Energy is left in the detector by the recoiling silicon nucleus. 
Thus the energy E^ must be found as a function of the neutron scattering 
angle v, and the form of the differential cross-section must be known, 
in order to compute the probability of an incident neutron leaving a 
countable energy in the detector. In particular, the values of v must be 
found for which .0015 MeV 5 E^ ^.010 MeV or equivalently, E^-.OIO 
MeV s =* Ej-. 0015 MeV (a low energy count) and .010 ^ E^ ^ . 030, 
or equivalently, E^-.03 £ E^ ^E^-.01 (a high energy count). The 
energies and angles are related by (Evans 8 ) 




Differential cross-section data are apparently not available 
for 14^ x j on the assumption that the difference between one element 
and its neighbor will not be very great, the data for ^Al were used. 7 



The results are summarized in Table III. The pulse height defect, 
due to nuclear (non-ionizing) collisions, will increase these numbers 
slightly . 8 

2 . Si^ (n,n l ) Si^^ (inelastic scatter) 

-2 

From Figure 11, the probability is <2 X 10 , so the total number 

2 2 

of interactions (using a total incident flux of 8.4 X 10 / cm sec) is 

2 28 
^17 cm sec. The first excited level in Si is 1.78 MeV above ground. 

Thus from Figure 2, the probability of leaving a count in the low energy 

• -4 -4 

detector is at most ^10 , and for the high energy detector, X 10 

The contributions to the counting rates are therefore 

(17)(1 0 -4 )(2) = . 0034 counts/sec, 1.5-10 keV 

(1 7)(2 x 10" 4 )(10) = .034 counts/sec, 10-30 keV 

3. Si (n, p) A1 ; A1 -* 0 + Si , T^ - 2.3 min 

_ 3 

From Figure 11, the probability of this reaction is <6 x 10 , 

and the threshold is 3. 86 MeV. For incident neutron energies in the 

range 3.86 ^ E ^ 3.96, the proton is left with E ^ 30 keV . 8 From 

Figure 9, a rough estimate of the number of incident neutrons in this 

2 

energy range is ^3 /cm sec. Thus as upper limits to the count rates 
from this reaction we have 

(3)(6 X 10 ^)(2) = . 036 counts/sec, 1.5-10 keV 
(3)(6 X 10~^)(10) = .18 counts/sec, 10-30 keV 

28 

The AT nucleus may be formed in an excited state; the first 

energy level is 31 keY above ground and up to excitation energies of 

7.7 MeV, de-excitation is by gamma emission . 9 The 31 keV level has 

-9 

a lifetime of about 2 X 10 sec; angular momenta for higher levels are 
uncertain, but de -excitation probably by dipole transition, so other 
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lifetimes are of this order of magnitude. Thus radiative de-excitation 
28 

of the A1 nucleus will be coincident with the energy loss of the proton, 

and no extra counts will result. 

28 

The beta emission of Al is followed by a radiative de-excitation 
28 

of the Si nucleus. These events may deposit counts in the detectors 
since they occur several minutes after the (n, p) reaction. 

I .*78 MfcV 

.s«io'“s 


Si 88 



The end point of the beta spectrum is 2.87 MeV ; thus a beta 
can deposit a count either by being ejected with an energy less than 
30 keV, or by being ejected with a higher energy near the surface of 
the detector, and losing less than 30 keV before it escapes. The beta 
spectrum should look like the one given by Evans 8 in his Figure 1.4, 
Thus a very rough estimate of the fraction of betas created with an 

_3 

energy less than 30 keV is X 10 . From Figure 9, the number of 

r 2 

incident neutrons above the threshold energy is approximately 30 cm 


sec; since the probability of a (n, p) reaction is about 6 X 10 ,. the 

number of betas created with E ^ 30 keV is (30)(6 X 10 ^)(2 X 10 = 

-4/2 

3.6 X 10 /cm sec. 



For betas with energies around 2 MeV, dE/dx — .4 keV /micron 10 
in silicon; thus a beta would have to traverse only 3-70 microns to lose 
1.5-30 keV. The probability of a (n, p) reaction in the top or bottom 
70 microns is (Figure 10) 

_ . 4 barns /. -6.4x10 ^ cm^/g xl.65 XlO ^g/cm^N 

3 barns \ e / . 

' 'SB. 2. 8 X 10" 4 
2 

Using 30/cm sec neutrons above threshold energy, the number of inter- 

_4 2 

actions producing a count is roughly 84 XlO /cm sec. 

An interaction of the 1.78 gamma, since it is coincident with 

the beta, would simply add to the energy deposited by the beta. Thus 

these events will be neglected. The total number of interactions producing 

-4/ 2 

counts is then 88 X 10 /cm sec and the contributions are 

.018 counts/sec, 1.5-10 keV 

.088 counts /sec, 10-30 keV 
28 

from the decay of the Al nucleus. Adding these to the (n, p) contri- 
butions on page 24, the totals are , 

.054 counts /sec, 1.5-10 keV 
.27 counts/sec, 10-30 keV 

4. Si (n, a) Mg I 

_ 3 

The probability of this reaction is ~ 10 (Figure 11), and the 

2 

threshold is 2. 66 MeV. Using an incident flux of 210/cm sec, 

/ 2 25 

21/ cm sec reactions take place. Mg is a stable nucleus; up to an 

energy of about 7 MeV above ground, de -excitation is by gamma emission, 

and the first excited state is .6 MeV above ground. 9 Lifetimes are of 
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_9 

the order of 10 seconds, so the gamma energy loss is coincident with 
the alpha. 

As above, the alpha can leave a countable energy either by being 

emitted with an energy of less than 30 keV, or by being emitted with a 

higher energy near the surface of the detector. 

For incident neutron energies 2.65 ^ E ^ 2. 9 MeV, the alpha is 

emitted with less than ~30 keV. From Figure 9, roughly 55/ cm sec 

neutrons are incident with these energies, so the contribution from this 

2 

type of interaction is . 055 /cm sec. 

At higher energies, dE/dx = 200 keV/micron for alphas; 10 to 
lose 30 keV, an alpha must traverse only 7 microns, and the probability 
of a (n,cy) reaction in the top or bottom 7 microns is approximately 

^ , . 1 barn /, -6. 4X10 2 cm 2 /g X 1. 65 X 10 3 * g/cm 2 \ 

2 * 3 barns V 1 ' e _ J 

27X10 " 5 

/ 2 

Using an incident flux of 210/cm sec above threshold, the contribution 

2 

is .015/cm sec. Thus the total number of counts/second is 
-(. 07) (2) = .14 counts/sec, 1.5-10 keV 
(. 07)(10) = .7 counts/sec, 10-30 keV 
Table IV summarizes the effects of neutrons in silicon. 

Table IV . .. 

Interaction Count Rates Due to Interaction 

, 1. 5-10 keV 10-30 keV 

(n, n) 7. 1 c/s 39 c/s 

(n,n‘) . 0034 c/s .034 c/s 

(n, p) .054 c/s .27 c/s 

(n,or) 


. 14 c/s 


7 c/s 
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C. RTG Neutrons: Indirect Effects Due to Lead Shield 

The following neutron reactions can occur in the lead shield: 


1. 

Pb 207 

207 

(n, n) Pb (elastic scatter) 

2. 

Pb 207 

207* 

(n, n 1 ) Pb (inelastic scatter) 

3. 

Pb 207 

, , 208 
(n,y) Pb 


The cross-sections for these reactions are shown in Figure 12, and 

2 

the probabilities for 3 g/cm of lead are shown in Figure 13. 

1. Pb 2 ^ 7 (n, n) Pb 2 ^ 7 (elastic scatter) 

This reaction has the effect of scattering neutrons, not headed 

for the detectors initially, into the detectors. Using the formula on 

page 23, and noting that the transparency, of the lead shield to neutrons 

2 

means that neutrons have two chances to scatter, . 4rt neutrons / cm 

2 

sec head for the detectors, where n is the number /cm sec which 

2 2 

elastically scatter. Since 8.4 X 10 / cm sec are incident, and the 

-2 2 

probability of an elastic scatter is ^4. 5 X 10 , 15/cm sec will be 

scattered into the detectors. Referring to Table III, the maximum 
contributions are: 

(1 5) (2) (. 73 X 10" 2 ) = .22 c/s, 1.5-10 keV 
(1 5) ( 1 0) (. 81 X 10‘ 2 ) = 1.2 c/s, 10-30 keV 

2. Pb 2 ^ 7 (n,n') Pb 2 ^^* (inelastic scatter) 

Since the probability of an inelastic scatter is at most about 
one -half the probability of an elastic scatter, the contributions from 
scattered neutrons will be about one-half of those given above: 

.11 c/s, 1.5-10 keV 
.6 c/s, 10-30 keV 
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207 

For excitation energies up to about 6.7 MeV, the Pb nucleus 
de -excites by gamma emission. From Figure 13, the number of 
inelastic scatters is <(8.4 X 10 )(2. 5 X 10 ) = 21 /cm sec. The first 
excited state is . 57 MeV above ground, so from Figure 2, the contri- 
butions are: 


(. 4)(21)(2)(8. 5 X 10 ) = .014 c/s, 1.5-10 keV 

(. 4)(21)(10)(2 X 10" 3 ) = .168 c/s, 10-30 keV 

Thus the total contributions due to neutrons and gammas from 
this interaction are: 

. 12 c/s, 1. 5-10 
.23 c/s, 10-30 


' on 207 / * -du 208 

3. PB (n, y) Pb 


The Q value for this reaction, using the mass defect values 

given by Leighton, 11 is 37. 79 + 8. 36 - 38. 77 = -7.4 MeV. De-excitation 
208 

of the stable Pb nucleus follows by gamma emission direct to ground 
in ^ 90% of the cases. 9 As there are evidently no neutrons with energies 
above 5 MeV, this reaction can be neglected. 

Table V summarizes the effects of RTG neutrons due to a lead 

shield. 

Table V 

Effects of RTG Neutron Interactions in a Lead Shield 

Counting Rates 

Interaction 1. 5-10 keV 10-30 keV 

(n,n) .22 c/s 1.2 c/s 

(n, n') . 12 c/s .23 c/s 
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D. RTG Neutrons: Indirect Effects in a Csl Shield 

1 27 

The following interactions can occur with the I nucleus: 

1. I* 27 (n, n) I 127 (elastic scatter) 

127 127* 

2. I (n, n') I (inelastic scatter) 

. ,127 , > t 128 t 128 - , v 128 ' , . . . 

3. I (n.,y) I ; I -• 6 + Xe (T^ = 25 minutes) 

2 

(Beta spectrum end point is 2. 12 MeV) 

The Iodine cross-sections are shown in Figure 14, reproduced from 

Hughes/ Cs cross-sections are evidently not available. Since the total 

2 

iodine cross-section is approximately 6 barns/atom = . 0285 cm /g, 
the shield is almost transparent to neutrons. Noting that this trans- 
parency gives most neutrons two chances to interact, and each passage 

2 2 
through the shield (2.9 g/cm thick) involves approximately 1.45 g/cm . 

_ p Q v 0? ft ft 

Iodine, the fraction which interacts is 1 - e — . 0825. Since 

the recoiling Iodine nucleus would cause an anticoincidence signal for 
most of the above reactions, interactions in the active shield would cause 
no extra counts in the detectors. However, the contribution to the dead 
time would be (,0825)(8.4 X 10 2 /cm 2 sec)(150 cm 2 ) = 10^ c/s. 

E. Diffuse Cosmic X-ray Leakage 

An upper limit to this leakage flux can be had by considering 

the total number of diffuse x-rays with energies greater than 200 keV 

since the shield becomes quite transparent at higher energies. From 

-12 

page 10, this is approximately 10 /cm sec sr. Using Figure 2, the 
number of interactions is 

$ (10’ 1 )(4tt)(2)(3 XlO' 3 ) = . 0024 c/s, 1.5-10 keV 
^(10' i )(4n)(10)(1.5 X10' 2 ) = .06 c/s, 10-30 keV 
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From page 10, the number of diffuse x-rays 10 <E < 200 is 3. 6/crn 

sec sr. Almost all of these would interact in an active shield, giving 

4 

a dead time contribution of 1.35 X 10 counts /second. 

F. Cosmic Ray Interactions in the Shield 

2 

Using a cross-section of .0125 cm /g for cosmic rays, a 

3 g/cm shield will cause a fraction (1-e ) = 7.5 X 10 

of incident cosmic rays to interact. Over 2tt sr, there are approxi- 

mately 6. 28 cosmic rays/cm sec and .47 interactions /cm sec. Using 

the argument of page 23 and noting that there are four walls of the 

2 

shield producing gamma rays, (. 8)(. 47) = . 38 gammas/cm sec are 

directed towards the detectors, if every cosmic ray interaction produces 

a gamma ray. Using the maximum probabilities of Figure 2 to obtain 

an upper limit, we have 

(. 38)(2)(6 X 10‘ 2 ) = 4.5 X 10‘ 2 c/s, 1.5-10 keV 

(. 3 8) ( 1 0) (3 X 10" 2 ) = .11 c/s, 10-30 keV 

2 2 

In an active shield, (,47)(600 cm ) = 2.8 X 10 interactions /sec would 
occur. 

G. Cosmic Ray Interactions in the Spacecraft 

2 

If N is the flux of cosmic rays (per cm sec sr), Nrra cosmic 
rays/sec will be incident on the area a bounding a volume V. 



If the volume V has a mass W, a typical path through the volume has a 

2 2 
thickness W/a g/cm . Using a mean free path of 80 g/cm for, .cosmic 

rays, the probability of an interaction is 1 - e Rough estimates 

I - i 
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5 4 2 

give W = 6.5 X 10 g and a = 6 X 10 cm for the Grand Tour spacecraft, 

2 . . 

so W/a = 11 g/cm , and the probability is approximately W/80a in this 
case. Thus the number of interactions per second is Nira W/80a = 
NrrW/80. Assuming uniform production of gammas throughout the volume, 
and assuming that each interaction produces a gamma ray, then a small 
volume dV produces y/sec. The flux at a detector a distance 

r from dV is — —■ y/ err? sec. A detector of area A, whose 

normal makes an angle 0 with r, will intercept 


NttW dV _1 (. 

80 V . 2 A V 

4irr 


NttW 

80 


Thus the number of counts /second due to the entire volume is 


NttW 1 A\ix ? dV 

80 V 4tt J 2 
r 

Evaluation of this integral for simple geometries (spheres and cylinders) 
shows that it is of the order 4rrr, where r is a characteristic dimension 
of the object. Taking the case of a sphere, with the detector at the 
center, there will be 

NttW 1 Ap,x _ NttW 
80 4 3 4 tt " 80 

3 " r 

2 

Using N = 1, W/a = 10 g/cm , this becomes 1.2 Apx counts/sec. 
Taking the maximum probabilities from Figure 2, we get 

(1. 2) (2)(6 XlO" 2 ) = .14 c/s, 1.5-10 keV 

(1. 2) (10)(3 X 10’ 2 ) = .36 c/s, 10-30 keV 


A - Ntt 
2 AM-x g0 


A p x 
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IV. SUMMARY 

The emphasis has been on the derivation of upper limits; a number 
of factors could reduce the counting rates. Among them are (1) the 
presence of matter between the RTGs and the detectors (for example, 
the propulsion bay might cast a "shadow" on the science end of the space- 
craft); and (2) the presence of a housing around the science portion of 
the spacecraft. While these might not do much for the neutron induced 
background, they would certainly reduce the gamma background. In 
addition, a housing might reduce the necessary shield thickness. 

To complete the comparison between the active and passive shields. 
Table VI summarizes the counting rates from all sources for an active 
and for a passive shield. To compute the factor B/(l-bd), a dead time 
gate length must be assumed. Under normal circumstances, a reasonable 
assumption might be that a 30 microsecond gate is generated whenever 
an interaction occurs in the shield. In the present case, however, the 
shield is counting a large number of interactions due to the RTGs, and 
only a small fraction of these interactions results in a count occurring 
in the main detectors. Thus it might be feasible to send the main 
detector pulses through a delay line, meanwhile analyzing the main 
detector and shield pulses (perhaps by means of a peak detector) for 
coincidence to within a microsecond or so, and generate a 30 micro- 
second gate only for true coincidence between the shield and the main 
detector. In this way, the amount of spurious dead time could be 

reduced to almost zero. As Table VI shows, the number of counts 

4 

per second in the shield is 4.1 X 10 , which is sufficient to cause 100% 
dead time with a 30 microsecond gate length. 
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Table VI 


« 

Passive Collimator 
(counts /second) 

Active Collimator 
(counts /second) 



Source 1 

.5-10 keV 

10-30 keV 

1.5-10 keV 

10-30 keV 

Dead Time 

Cosmic X-rays 
(Forward Aperture) 

.53 

.55 

. 52 

. 52 

- 

- 

RTG Gammas 

2.3 

12 

1.6 

9 

1.8 

4 

X 10 

RTG Neutrons 
■; (Direct Effects) 

7.3 

40 

7.3 

40 

__ 


RTG Neutrons 
(Indirect Effects) 

.34 

1.5 

0 

0 

10 4 


Cosmic X-ray 
Leakage 

. 0024 

.06 

~ 0 

/~0 

1.3 

o 

f— < 

X 

Cosmic Ray 
Interactions 
in Shield 

. 045 

.11 

0 

0 

2. 8 

X 10 2 

Cosmic Ray 
Interactions in 
Spacecraft 

.14 

.36 

14 

36 

~ 0 


TOTALS 

10.65 

54. 58 

9.56 

49. 88 

4. 1 

x io 4 

TOTALS WITHOUT 

NEUTRON 

CONTRIBUTIONS 

3. 01 

13. 08 

2.26 

9.88 
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The total background counting rates are computed with all inter 
actions included, and with all interactions except neutron interactions 
included; since neutron interactions are responsible for most of the 
background counts, some sort of pulse shape discrimination should be 
investigated to reject this contribution. 
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PART III 


DESIGN STUDY FOR AN X-RAY EXPERIMENT 
USING A PROPORTIONAL COUNTER 
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I. THE DETECTOR 

The primary purpose of the Grand Tour X-ray Experiment 
described in this report is to measure the spectrum and time variations 
of auroral zone x-rays in the range 2-30 keV, and the size and location 
of the auroral zones of Jupiter and other outer planets. The secondary 
purpose is to measure the spectra of cosmic x-ray sources and their 
time variations. , 

This design study is based on the Reuter-Stokes proportional 
counter RSG-77 (Figure 1) filled to 1 atmosphere with 97% Xenon, 3% 
CO^. and with a 5 mil Beryllium window. This detector must be 
collimated, both to avoid the problem of source confusion while viewing 
an x-ray source, and to provide a characteristic response pattern which 
will allow the angular size of an x-ray source to be measured. The 
first requirement is satisfied by rough collimation to a field approxi- 
mately 5° X 5°. In order to measure the size of Jupiter's auroral 
zones before the spacecraft enters the magnetosphere (and the detector 
background counting rate increases drastically), much finer collimation 
is needed: since the Jovian magnetosphere probably extends out to 

50 Rj in the solar direction, Jupiter will subtend only about 1° when 
the spacecraft enters the magnetosphere, and collimation down to 
ro 1/10° Will be necessary. This suggests the use of two collimators; 
a simple slat collimator for the 5° X 5° field, and a wire grid modulation 
collimator for the measurement of angular sizes (Figure 2). Consider 
first the slat collimator, which divides the active area of the proportional 
counter into an array of squares. If there are n slats/unit length, 
and their height is h, then each square is collimated to an angle 
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8 = tan ^(l/2nh). Since there are (n£+l) slats along the length each 

having dimensions w X h, and nw+1 slats along the width, each having 

dimensions 1 X h, the total area of the slats is (ni + l)wh + (nw + 1) 

fh = 2n f-wh + h(w + 1) = + (' 2 I tanv) > Thus » for a g iven pro- 

portional counter (i, w) and collimation angle 8, the minimum total area 
of slats is obtained for the maximum n. A minimum total area is 
desirable not only because it minimizes the weight, but also because 
it keeps the mass in front of the active area down, thus minimizing 
scattering into the proportional counter. 

For the RSG-77, a natural choice for n is 1/1.14 cm since it 

allows the collimator slats to be attached to the Beryllium window 

support structure. Thus for 8 = 5°, h = (2n tan 9) * = 5 cm, and the 

2 

total area is A = 757 cm . 

Referring to the modulation collimator response pattern shown in 
Figure 2, we see that if the characteristic angle cp = 1/5°, a distant 
point source transiting the collimator will cause the proportional counter 1 
to go from a maximum to a minimum count rate in 1/10°. It follows 
that, when the x-ray source begins to deviate from a point source and 
subtends an angle of 1/10°, the response pattern will start to change 
noticeably. Figure 3 shows the efficiency of tungsten wire for stopping 
x-rays; choosing 8 mil wire, and noting from Figure 2 that the wire 
spacing should be approximately equal to the wire diameter in order 
to block the x-rays coming from angles |>cp, -|cp, etc., the height of 
the collimator b = a/tanco = 6. 1 cm. Although the modulation collimator 
is shown in Figure 2 as a separate unit, it could be incorporated into 
the slat collimator to save space and weight. 
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It will be assumed that all interactions in the proportional 
counter which are not in the 2-30 keV energy range, as well as 
interactions in the slat collimator, if it is active, will give rise to 
an anticoincidence signal. (An active "slat" collimator might consist, 
for example, of a Cesium Iodide slab, viewed by PM tubes, with an 
array of cylindrical holes leading to the proportional counter.) In 
addition, interactions in the proportional counter which are in the 2-30 
keV window, but are not caused by the x-ray source being viewed, 
will contribute to the background counting rate. The relative importance 
of these two effects is as follows. Suppose that an x-ray source, 
whose true strength is S counts/sec, must be detected in the presence 
of a background count rate whose true strength is B counts/sec. The 
detector dead time is bt, where b is the number of interactions/sec 
which produce an anticoincidence pulse, and t is the length of the gate 
signal for such a pulse. Then, in a time interval T, ST-STbt counts 
are collected from the source, and BT-BTbt from the background. 

The source may be considered to be detected if ST(l-bt) s 3a, 


9B 


bkgd 

, and the factor B/(l-bt) determines 


3 /BT(l - bt) . Thus T ^ §^"(1 - bt) 
the time needed to detect a given source at a 99% confidence level. 

Figure 4 shows the x-ray absorption coefficient, p,, for Xenon, 
versus x-ray energy. 1 ’ 2 It is important to note that, up to an energy 
of approximately 300 keV, the photoelectric absorption coefficient is 
much greater than the Compton coefficient; this means that up to 300 
keV, if an x-ray interacts at all in the proportional counter, it will 
most probably interact photoelectrically, leaving all of its energy. 

Thus x-rays in the energy range 30-300 keV which are incident upon 
the detector will contribute to the anticoincidence count rate b, rather 
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than the background count rate B. 

This is shown in greater detail in Figure 5, which gives, for 
the RSG-77 counter, the probabilities of various types of interactions. 

The top graph (A) gives the probability of either a photoelectric or a 
Compton interaction; that is (probability of no interaction in . 009' 
Beryllium) X (probability of interaction in 1" Xenon at STP) = 
exp(-p n x„ ) [1 - exp(-(j, x )], where and u are the attenuation 

JjG dG X6 X6 IjG X© 

-3 2 

coefficients for Beryllium and Xenon, and x^ e = . 005" = 9 X 10 g/cm , 

-2 2 

x xe = 1" = 1.37 X 10 g/cm . The lower graph gives the probability 
that a photon, if it interacts, will leave between 2 and 30 keV in the 
detector. In the energy range 2-30 keV, this is the same as the top 
graph. Above 30 keV, this is 


(Probability of interaction of any kind) X 


Cross-section for transferring 

2-30 keV to electrons 

Total cross-section 


1 ^ 

= e X p<-^ x ) [1 - expf-^ X )] — f 

total 2 


where 0^^ is the total cross-section for a photon interaction of any 
kind, and ■ ^ Compton ^^ eren ^ a ^ Compton cross-section for 
a Compton interaction which leaves the electron with an energy 
T -• T + dT; this integral was done numerically, using the graphs of 
Nelms. a Thus the lower graph gives the probability that an incoming 
photon contributes to the background rate B, while the difference 
between the top graph and the lower graph gives the probability that 
an incoming photon contributes to the anticoincidence rate b; this 
latter probability is zero for energies 2-30 keV and to a good approxi- 
mation, equal to the top graph for other energies. 
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In the following sections the counting rates b and B will be 
determined for this detector; it is convenient to distinguish between two 
different types of interactions. 

1. Counting rates through the forward aperture, caused by 
particles and photons which enter the detector after 
passing through the two collimators. 

2. Leakage counting rates, caused by photons and particles 

which originate outside the angles of collimation, but 
"leak" through the collimator or the sides of the pro- 
portional counter. > 

II. BACKGROUND COUNTING RATES THROUGH THE FORWARD 

APERTURE 

A. Diffuse Cosmic X-ray Flux 

The spectrum of the diffuse flux is 4 

jg = 10 E **^/cm^ sec sr keV, 1 <E <10 keV 
= 16. 8 E ** ^ /cm^ sec sr keV, 10 <E < 40 keV 
= 128 E cm^ sec sr keV, 40 < E < 1000 keV 

This spectrum is plotted in Figure 6, along with 

1. the product of the spectrum and the probability of either a 
photoelectric or Compton interaction (top graph of Figure 5), 
and 

2. the product of the spectrum with the probability of an 
interaction which leaves 2-30 keV (lower graph of 
Figure 5). 
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Integration of these two graphs gives 

2 

1. 7.15/cm sec sr, the flux which produces the 2-30 keV 

background counting rate, and 
2 

2. . 162/cm sec sr, the flux which produces the >30 keY 
anticoincidence counting rate. 

2 

Using a total effective area of 39.2 cm for the RSG-77, and a 

field of view of 5° X 5° = 7. 6 X 10 ^ sr, and noting that the modulation 

collimator obscures half of the counter's effective area, the geometrical 

2 

factor relating flux to count rate becomes . 149 cm sr for the 2-30 
keV background count rate. (Although the modulation collimator 
becomes transparent at energies above 30 keV, the contribution to the 
2-30 keV count rate from fluxes over 30 keV is negligible. ) Thus 
B = 1.06 counts/sec. 

Since the anticoincidence count rate is caused by photons with 

energies greater than 30 keV, an upper limit to this counting rate may 

be made by considering the modulation collimator to be transparent; 

2 

the geometrical factor is therefore . 398 cm sr, and the count rate 
-2 

is b = 6.45 X 10 counts/sec. 

HI. LEAKAGE FLUXES 

A. RTG Gamma Rays 

The spectrum of RTG gammas is shown in Figure 7; this is 
the spectrum at 3 ft, 0° from the RTG given by Noon et al .f 
multiplied by 3 and extrapolated to 13 ft, to give the effect of 3 RTGs 
assuming no matter exists between the RTG and the x-ray experiment. 
To obtain the counting rates caused by this flux, consider first the 
effect of the orientation of the detector to the RTG (Figure 8). To 
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a first approximation, there can be three orientations: 



Figure 8 

1. The side H XL faces the RTG; photon path length through 
counter is W. Count rate « HL [1 - exp(-p W)] . 

2. The side H X W faces the RTG; photon path length L. Count 
rate « HW [1 - exp(-p, L)]. 

3. The side L X W faces the RTG; photon path length H. Count 
rate « LW [l - exp(-jj, H)]. 

Noting that the RTG gives no flux below 50 keV, and referring to 

2 

Figure 4, the absorption coefficient for Xenon is 12 g/cm at 50 keV. 

~Ui X 

Using the approximation 1 - e ^ — px, it follows that all of the orien- 

tations should give approximately the same count rate. As a worst 

-2 2 

case, the length L, = 6.25 inches = 8 . 55 X 10 g/cm , so that, at 
50 keV 1 - e = 1 - e = ^ 34 ^ w hile pL = 1.02. However, 

—X • . 

note that since e ^ 1 - x, the use of this approximation will give 
an upper limit to the counting rate even when it is not a good approxi- 
mation. 

Also plotted on Figure 7 are the products of the RTG flux with 
the maximum probabilities of Figure 5 for each energy range; thus the 
AA curve gives the flux responsible for the interactions which leave 
>30 keV with an electron, while the +++ curve gives the flux responsible 
for the 2-30 keV interactions. Integrating these curves, and using the 
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2 

total active area of the detector, 60.5 cm (since the aluminum support 
structure is transparent to x-rays), the counting rates are: 

B = 24 counts/sec 

2 

b = 5.3 X 10 counts/sec 

The previous analysis neglects the following possibility. 

An electron with energy greater than 30 keV, generated either 
by a Compton or photoelectric interaction, could leave 2-30 keV in 
the detector, if its path through the gas were less than its range. 

This is illustrated in Figure 9» which plots the distance which an 
electron must travel to deposit 2 and 30 keV, versus electron energy 
for Xenon at STP. Thus, for example, a 1000 keV electron, traveling 
1 cm through the detector would leave a countable energy. On the 
other hand, it is evident from Figure 7 that the worst offenders by 
far are the RTG photons in the 100-200 keV energy range, and it 
will be shown later that these can be shielded against relatively easily. 

RTG photons can produce background counts in yet another way, 
illustrated in Figure 10: an incoming photon can interact in the 

counter walls, either photoelectrically or Compton scattering, thus 
producing an electron which can enter the counter gas. 


-A |— t 



CROSS SECTION Of PROPORTIONAL CENTER 
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If the incoming photon has an energy E, it can produce an electron 
whose maximum energy is (E - binding energy of atomic electron) — E, 
since electron binding energies are small compared to RTG photon 
energies for materials like Aluminum and Beryllium. Let the range 
of the electron be R(E) in the counter wall material, and let the wall 
thickness be t. Then, if the electron is to enter the gas before losing 
all its energy in the wall, it must be produced a distance R(E) or less 
from the gas. Thus the incoming photon must not interact in the wall 
for a distance t - R(E) [probability exp-p (E)(t - R(E))] , and then inter- 
act in a distance R(E) [probability 1 - exp(-p(E) R(E))]. If it is 
assumed that all such interactions produce a count in the detector, then, 
as an upper limit the counting rate will be proportional to 

A e -p(E)[t-R(E)] j-j _ e -n(E) R(E)j 

where A is the cross-section area presented to the RTG photon flux. 

A comparison of references 2 and 6 shows that for all types of 

materials the product p(E) R(E) <<1 at RTG photon energies; thus the 

counting rate takes the form Ae p(E) R(E). The graph of 

Figure 11 gives p(E) R(E) e versus energy, for 1.15 g/cm^ of 

2 

Beryllium (.25" thick) and 1.15 g/cm of Lead (.04" thick). Thus, 

with the total wall weight held constant, the effect of high Z and low Z 

materials may be seen on this graph. Consider the Beryllium curve 

first. At energies over 50 keV, the photoelectric cross-section is 

negligible, and the product p(E)t is small enough so that the exponential 

term is close to one; thus the expression simply reduces to p(E) R(E) 

2 

(independent of wall thickness). The range R(E) « E , approximately. 
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while the Compton cross-section decreases slowly with energy; the 

2 

result is a curve which increases roughly as E , turning over at 
higher energies. 

The case for lead is considerably different; at low energies, 
the photoelectric cross-section is very large, so the exponential term 
is small, and the product is small up to approximately 100 keV. 

Then in the 100-200 keV region the photoelectric cross-section decreases 
becoming comparable to the Compton cross-section, bringing the expo- 
nential term close to 1, and allowing the p(E) R(E) term to dominate 
the expression; since p for lead is still larger than p for Beryllium at 
these energies, the total probability is greater for lead than Beryllium. 
As the energy increases further, p for lead becomes about the same 
as p for Beryllium (which simply indicates that the Compton cross- 
section dominates, and there are approximately the same number of 
electrons /gram for all materials) and the two curves are approximately 
equal. 

Since it is easy to shield against incoming RTG photons with 
energies below 200 keV, the use of a high Z wall material is not 
necessary, and Beryllium appears to be the best choice. 

RTG photons can produce background counts in yet another 
manner, illustrated in Figure 12: a photon can pass through the 

first wall and the gas of the counter without interacting, then back- 
scatter in the second wall and interact in the gas. A rough upper 
limit to the probability for this type of event is: (Probability that photon 
of energy E does not interact in the first wall)(Probability that photon 
of energy E does not interact in gas)(Probability that photon of energy 
E Compton scatters in second wall) (Probability that 180° backscattered 



49 . 

photon of energy E' interacts in gas). This product is shown in the 

2 

graph of Figure 13, for 1.15 g/cm Beryllium, and a 1.5" path 
through Xenon at STP. 


FI6. M 


Finally, the graph of Figure 14 shows the RTG gamma fluxes 

at the detector (same as Figure 7) and the product of these fluxes 

with the probabilities of Figure 11 (Beryllium curve) and Figure 13. 

2 

Multiplying the total of the two lower curves (17 x-rays/cm sec) 

2 

by the cross-sectional area (on the order of 40 cm ) gives an upper 

limit to b + B of 680 counts/sec. Thus these effects are potentially 

as serious a problem as direct interactions of the RTG photons. 

It seems likely that the propulsion bay of the spacecraft will 

shield the scan platform from the RTGs ; the amount of material present 

2 

is difficult to estimate, but, as an order of magnitude, 40 g/cm might 
be used (this corresponds to 15 cm of Aluminum). The effect of 
Compton scattering of photons is illustrated in the graph of Figure 15; 
from reference 3, graphs of Compton cross-section versus photon 
energy were combined to give the average energy of a Compton 
scattered photon versus distance traveled. At higher energies, where 
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the differential cross-section peaks in the forward direction, the 
distance traveled is more or less in a straight line; at lower energies 
the cross-section becomes roughly isotropic, and the distance traveled 
by the photon is not simply related to the thickness of the shielding 
material. However, some important features can be deduced from 
this graph: 

1. The 1-2 MeV RTG photons, which are potentially an impor- 
tant source of background, will be reduced in energy to 
roughly 200-450 keV. 

2. The 100-200 keV photons will be reduced in energy to 
roughly 50-70 keV. 

Comparing the photoelectric and Compton cross-sections in Grodstein 2 
it is apparent that Aluminum (in which the photoelectric and Compton 
cross-sections become equal at 50 keV) will do little to remove the 
RTG photons, although it will alter the spectrum significantly. To 
reduce the detector background, a material like Lead will be neces- 
sary (where the Compton and photoelectric cross-sections become 
equal at about 500 keV) or an active collimator. 


The percentage of photons which photoelectrically interact in a 

thickness x g/cm 2 of Lead is [pp E (E)]/[p, total (E)](l - e total ) 

where ypj,(E), are the photoelectric and total attenuation 

coefficients in Lead. From Grodstein 2 it is clear that for energies 

below 200 keV this percentage is greater than 75%, for 1 cm (11.34 
2 

g/cm ) of Lead, and reduces to about 40% at 500 keV. It seems likely 
that, after the backscatter calculations have been refined past the 


upper limit stage, 1 cm of lead will provide the necessary shielding 

2 3 

against RTG photons (1 cm X 40 cm x 11.34 g/cm — 450 gm). 
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This may be compared to the weight of an active shield as 
follows: assume that any interaction in an active shield, either photo- 

electric or Compton, gives rise to an anticoincidence signal. Then 

~^total X 

the attenuation factor for photons is 1 - e .. If we require the 

- ^total X 

same attenuation at 500 keV as for the Lead shield, 1 - e = .40 

2 

and g,x = .51. For Csl, this would give x = 7.6 g/cm , so the 

2 

weight would be 7.6 X 40 cm ss 300 gm. Considering the additional 
electronics needed for an active shield, it appears that the active and 
passive shields would have comparable weights. The actual decision 
between the two will have to await more detailed calculations; it should 
simply be noted here that shields of moderate weight are capable of 
reducing the background to an acceptable level. 

One final type of RTG gamma interference problem should be 
considered here: gammas which scatter in the slat collimator and 

are redirected towards the proportional counter. The geometry for 
this type of interaction is shown in Figure 16. 



K ^-Active AeeA 
dUTAWAV 6>iPE VIEW 

Figure 16 
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Assume that the scattered photons are emitted isotropically, and let 

2 


a.-ii 

! s 


N be the number which escape /cm sec. Then a thin strip of area 
dx/n emits Ndx/n photons/sec, and the flux at a distance J^x^ + ("^) J 
(the distance to the center of the proportional counter cell) is 
approximately 


Ndx 




The proportional counter cell, of area A presents an effective area 

A coscp to this flux, and the probability that a photon will interact in 

cosC P 

the detector is 1 - e , where d is the detector depth. 

Thus, as an upper limit to the contribution from the thin strip, we have 


A Ndx 

. f 2 /I \2“i ,fJ Xe 

4 nn i x + fer) J 

A N ^X e d -1 

Integrating over x from 0 to h, this becomes tan 2 nh. 

2 

Using A = 1.308 cm , and 2 nh = 8.75, this reduces to .3 N ^ e d/sec. 

Consider next the effect of the orientation of the collimator. 
Figure 17 shows the detector after rotation about an axis normal to 
the Be windows; two sides make an angle 9 with the incident flux, and 
two sides make an angle 90-0; thus if each side has an area a, the 
number of interactions in the sides is proportional to 

, -^ol t/c ° se ') . 

2a cos 9 \1 - e / ■ + 

f -p r ,t/cos(90-9)\ 

2a cos(90-9) Vl - e / ^ 4a ^ o jt 
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TOP VIEW 


Figure 17 


where is the attenuation coefficient for the collimator material. 

2 

Thus if the incident RTG flux is F(E) photons /cm sec in a given 


energy interval, N(E) = 4F{E) p,^ o ^(E)t. Since there are 30 cells, 
the detector counting rate is 36 F(E) p ( -. o j(E)t |j^ e (E r )d,- where E' is 
the energy of the scattered photon. Using the Xenon interaction 
probabilities of Figure 5, and choosing E' to be the energy of the 
85° scattered photon (i. e . , the minimum energy photon, to give the 


maximum interaction probability), the results are shown in Figure 18. 

2 

A .02 cm thick (. 227 g/cm ) Lead collimator has been assumed. 

The total counting rates are 

2 

b = 3. 03 X 10 counts/sec 
B =1.22 counts /sec. 

Table I summarizes the effects of RTG gamma rays. 


Table I 

Direct Wall Collimator 

Interactions of Interactions of Interactions of 

Source Diffuse X-rays RTG Gammas RTG Gammas RTG Gammas 

B(2-30) 1. 06 c/s 24 c/s 1.22 

— t <6. 8 x 10 2 c/s 

6.45xl0 _2 c/s 5.3 X 1 0 2 c/s 3.03 X 10 2 


b(>30) 
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B. RTG Neutrons: Direct Effects in Proportional Counter 

The RTG neutron spectrum is shown in Figure 19; this is the 

spectrum given by Noon et al . 6 , multiplied by three, and extrapolated 

to 13 ft, as for the gamma spectrum. The total flux of neutrons is 

1.37 X 10 /cm sec. From Hughes, 7 it is clear that total neutron 

cross-sections are relatively constant above 10 keV and do not differ 

considerably from one element to its neighbor; from Walt and 

Barschall 8 , the total cross-section for Xenon can be estimated at 

-2 2 

6.5 barns/atom = 3 X 10 cm /g. By the same argument as in 
section HIA, the detector orientation does not affect the number of 
interactions, to a good approximation. Thus, considering the neutrons 
to be incident on the window of the proportional counter, the number 
of interactions /sec is 

, „ u , a 2 / 2 2 -3 x 10' 2 X 10' 2 . , , , 

1.37 X 10 / cm sec X 58 cm X(l-e ) = 3. 3/sec. 

The neutron will leave energy with the recoiling Xenon nucleus. From 
Evans y (p. 412), the maximum energy that can be imparted to a Xenon 
nucleus by a 7 MeV neutron is about 200 keV. The Xenon nucleus 
then loses energy both by ionization and nuclear collisions (Rutherford 
scattering by a screened potential). As the energy lost in the latter 
process involves no ionization, it is not detected in the proportional 
counter. From the pulse height defect curves given by Bertolini, 10 
it can be estimated that only Xenon nuclei with energies between 100 
and 200 keV will leave a countable energy; neutrons capable of 
transferring this energy must have an energy greater than about 
3 MeV, and the total number of neutrons above this energy is about 



NEuTgoMS/ CM 2 SEC 
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2 

13/cm sec, or 1/10 the total flux. Thus a rough estimate of the 
number of counts produced by direct neutron interactions is . 33/sec. 


C. RTG Neutrons: Interactions in the Proportional Counter 

Body 

Neutrons interacting in the proportional counter body can produce 

heavy ions and nuclear gamma rays which can go on to interact in 

the detector. Suppose, as in section IHA, that the counter body is 
2 

1.15 g/cm Beryllium. The Beryllium cross-section for neutrons is 

shown in Figure 20. Consider first the background due to Beryllium 

ions; these will have a range which can be estimated at 10 mg/cm , 1X 

for 10 MeV ions in Beryllium. Thus an upper limit to the probability 

of an ion being created and reaching the active volume of the counter 

2 

is (probability that incoming neutron goes (1.15-.01) g/cm without 

2 

interacting$(probability that neutron interacts in .01 g/cm ). Taking 

2 

the cross-section as 6 barns = .405 cm /g, this is 


e"l*14 X ,405 (1 - e - ' 405 x • 01 ) = . 63 X 4. 05 X 10 -3 = 2. 55 X 10 -3 


2 2 

Using the total flux of 1.37 X 10 neutrons/cm sec incident on a 

2 

side of the counter 40 cm in area, the total contribution is 14/sec; 

this represents an upper limit to the sum b + B. 

The background due to gammas can be estimated by noting that 

the first excited level of the Be nucleus is about 1.75 MeV 12 above 

ground. Thus the probability of a gamma being created and leaving 

a count in the detector is (probability of neutron interaction in 1.15 
2 

g/cm Be)(probability that 1.75 MeV gamma interacts in detector). 
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2 

Again, for the total flux of neutrons incident on a 40 cm side of the 
detector, an upper limit is 

b = (1 - e . A )(5 X 10 ) 1. 37 X 10 X 40 = (. 628)(5 x 10 *) 

X 1.37 X 10 2 X 40 = 1.7 c/s 

-40S y 1 1 ^ -S ? 

B = (1 r- e )(1 0 )(1 . 37 X 10 )X 40 = .036 c/s 


D. RTG Neutrons: Effects Due to Interactions in Lead 

Shields and Collimators 

The following neutron reactions can take place in lead for neutrons 
with energies up to about 7 Me Vi, 

1. Pb 2 ^ {n, n) Pb 2 ^ (elastic scatter) 

2. Pb 207 (n,n‘) Pb 2 ^^" (inelastic scatter) 

The cross-sections for these reactions are shown in Figure 21. 7 

In the slat collimator, an elastic scatter could direct a neutron 
towards the proportional counter, while an inelastic scatter could pro- 
duce a gamma (minimum energy approximately 570 keV 13 ) which could 
interact in the proportional counter. By the arguments of section HIA, 

upper limits to the counting rates produced by these effects are 

2 

36 F(E) (Probability of neutron interaction in . 227 g/cm Pb) 

_2 

(Probability of neutron interaction in 1.37 X 10 
2 

g/cm Xe) for neutron elastic scattering, and 

2 

36 F(E) (Probability of neutron interaction in . 227 g/cm Pb) 

_ 2 

(Probability of 570 keV gamma interaction in 1.37 X 10 
2 

g/cm Xe), for neutron inelastic scattering, where F(E) 


refers here to the total flux of neutrons. 


i 
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-2 2 

Using 10 barns = 2.91 X 10 cm /g for the cross-section of lead, 

-2 2 

and 6.5 barns = 3 X 10 cm /g as the cross-section for Xenon, the 
first count rate is 

36(1. 37 X 10 2 )(1 - e“ ,22? X 2,91 X 1° )(1 - e' 1 * 37 X 10 2><3 X 10 

= . 0134 c/s . 


And referring to Figure 5 the second is at worst 

36(1. 37 X 1 0 2 )(1 - e‘* 227 x2 * 91 x 10 )(1. 3 X 10" 3 ) 

If a 1 cm thick lead shield were placed along one side of the 

counter, as suggested in TIT A , for shielding against gamma radiation, 

neutrons which scattered inelastically in the lead could give rise to 

gammas which could then interact in the counter. An upper limit to 

the probability of this type of interaction is (probability that a neutron 

interacts in 1 cm = 11.34 g/cm lead) (probability that a 570 keV 

gamma interacts in proportional counter). Using the total cross-section 

2 

for lead, and the total neutron flux incident on 40 cm , we have 


= 4.24 X10' 2 


b = 40 X 1. 37 X 10 2 (1 
= 5.2 c/s 

B = 40 X 1. 37 X 10 2 (1 
= . 296 c/s 


e -U. 34 X 2 . 91 X10- 2 )(K5 x10 - 3) 

e -n.34 x 2 - 91 xl0 ‘ 2 )(7.5 X10- 5 ) 


E. RTG Neutrons: Effects in a Csl Shield 

2 

If an active shield consisting of 7.6 g/cm of Csl is used instead 
of a lead shield, neutron interactions in the shield will contribute to the 
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dead time counting rate b. From Hughes, 7 the cross-section of Iodine 

can be taken as 6. 5 barns /atom and from Walt and Barschall, 8 the 

cross-section of Cesium can be estimated at 6. 5 barns /atom. The 

-2 

attenuation coefficient for neutrons in Csl then becomes 1.5 X 10 

2 2 
cm /g, and the number of interactions/sec in a 40 cm shield is 


b = 40 X 1. 37 X 10 2 (1 - e" 7, 6 X lm 5 X 10 ) = 5.75 X 10 2 c/s . 

1 1 

Table II summarizes the effects of RTG neutrons on the proportional 
counter. 

Table II 


Source 

l l 


B 

Direct Interaction 
of RTG Neutron in 
Proportional Counter 

0 


. 33 c/s 

(n, n) Interaction 
in PC Body 

14/sec 

(n,y) Interaction 
in PC Body 

1.7 c/s 


. 036 c/s 

Neutron Interaction 
in Lead Collimator 

5. 58 X 10“ 2 
. c/s 

Neutron Interaction 
in Lead Shield 

5.2 c/s 


.296 c/s 

Neutron Interaction 
in Csl Shield 

0 


5.75 X 10 2 
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F. Diffuse Cosmic X-ray Flux Leakage 

The photons of the diffuse cosmic x-ray flux can "leak" through 

the slat collimator at angles greater than 5° (it is assumed here that 

the sides and bottom of the proportional counter will be surrounded by 
2 

several g/cm of electronics and mechanical structures, minimizing 

leakage through these areas). Assuming that the slat collimator is 

2 

made of .227 g/cm lead (section IHA), and referring to Grodstein 2 

it is clear that this collimator is at least 90% efficient up to 40 keV. 

Integrating the curves of Figure 6 above 40 keV, and assuming that 

2 

the entire top of the counter (~40 cm ) is exposed to 2tt sr, (geometrical 
2 

factor 40rr cm sr), the upper limits 

2 -1 2 

b = 60 cm X rr sr X 10 /cm sec sr = 18/sec 
2 -3 2 

B ■ = 60 cm X tt sr X 10 /cm sec sr = .18/sec 

G. Cosmic Ray Interactions in the Proportional Counter 

A minimum ionizing proton (energy greater than about 1 GeV) 

2 

will lose roughly 2 MeV/g/cm in Xenon. Paths from 0 to 7 inches 

are possible, so energy losses of 0 to 200 keV can occur. Using 

2 

a flux of 1 cosmic ray/cm sec sr above 1 GeV, and considering the 

2 

entire area of the counter (200 cm ) to be exposed to 2rr, the number 
of cosmic rays traversing the counter is 628/sec. Rise time discri- 
mination must therefore be used to minimize the effects of cosmic 


rays. 
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H. Cosmic Ray Interactions in the Spacecraft 

-2 

If N is the flux of cosmic rays (cm sec sr), N rr a cosmic 
rays /sec will be incident on the area a bounding a volume V. 



If the volume V has a mass W, a typical path through the volume 

!0 g/crr\ 
-W /80a 


2 2 
has a thickness W /a g/cm . Using a mean free path of 80 g/cm 


for cosmic rays, the probability of an interaction is 1 - e 

5 4 2 

Rough estimates give W = 6. 5 X 10 g and a = 6 X 10 cm for the 

2 

Grand Tour spacecraft, so W /a = 11 g/cm , and the probability is 
approximately W/80a in this case. Thus the number of interactions 
per second is Nna W/80a = NrrW/80. Assuming uniform production 
of gammas throughout the volume, and assuming that each interaction 
produces a gamma ray, then a small volume dV produces 

Y/ sec * The flux at a detector a distance r from dV is 

~ ~ a — ~~ p— v/cm^ sec. A detector of area A, whose normal 
80 V 4 tt r 2 ' 

makes an angle 9 with r, will intercept 


dV 1 AQ _ -px/cose. 0S NnW dV A p x . . 

80 V 4tt r 2 ’ 80 V 4n r a V /Sec * 
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Thus the number of counts/sec due to the entire volume is 


Nn W 1_ A u x .• dV 

80 V 4 tt Jr 2 


Evaluation of this integral for simple geometries (spheres and cylinders) 
shows that it is of the Order 4tr r, where r is a characteristic dimen- 
sion of the object. Taking the case of a sphere, with the detector at 
the center, there will be 


NnW 

80 



A u x 

4 TT 


4 tt r = 


3Ntt 

80 



2 

Using N = 1, W/a = 11 g/cm , this becomes 1.2 A p, x counts/sec. 

2 

Taking A = 40 cm and px = T for an upper limit, we get 48 counts/ 
sec for an upper limit to b + B. 
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I. Summary 

The following table summarizes the background effects from 
various sources. 


Table IV 


-• 

Source 

b 

B 


Diffuse x-rays 
(HA) 

.064/sec . 

1. 06/sec 


Direct Interactions 
of RTG Gammas 
(IIIA) 

530/ sec 

24/sec 


Interaction of RTG 
Gammas in PC 
Walls (IIIA) 

680/sec 



Interaction of RTG 
Gammas in 
Collimator (IIIA) 

303 c/s 

1 . 22 c/s 


Direct Interaction 
of RTG Neutrons 
in PC (HIB) 

0 

. 33 c/s 


(n, n) Interaction in 
PC Body (Hie) 

14 c/s 



(n,y) Interaction in 
PC body (IHC) 

1.7 c/s 

. 036 c/.s 


Neutron Interaction 
in lead Collimator 
(HID) 

. 055 c/s 



Neutron Interaction 
in lead Shield (HID) 

5. 2 c/s 

.296 c/s. 


Neutron Interaction 
in Csl Shield (HIE) 

575 c/s 

0 


Diffuse X-ray 
Leakage (HIF) 

18 c/s 

.18 c/s 


Direct Interaction 
of Cosmic Rays (HIG) 

628 c/s 



Cosmic Ray Interaction 
in Spa cecr aft (HIH) 

48 c/s 
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IV. CONCLUSION 

The most serious background problem is caused by RTG gamma 
radiation. Direct interactions of RTG gammas can be reduced by 
active or passive shielding of the counter, as well as pulse risetime 
analysis on the counter output; this technique can reduce the background 
by up to 90%, while retaining 80% of the x-ray efficiency. 13,14 

The problem of gamma interactions in the proportional counter 
body can be solved by the use of multiwire proportional counters: an 
array of anode wires is placed inside the counter next to the walls, 
and pulses are accepted from the central anode wire in anticoincidence 
with the wire array. Thus electrons originating in the counter walls 
are gated out. This technique has been investigated by Charpalc et al. 1 b 
Borkowski and Kopp, 16 and Bouclier et al . 17 In particular, Charpak 
et al . have shown that wire spacings down to 1 mm are possible with 
each wire behaving as a separate proportional counter, and with no 
sacrifice in energy resolution. Multiwire counters will also be used 
on the HEAO spacecraft. 

We conclude that, using currently produced proportional counters, 
in conjunction with proven anticoincidence techniques, an x-ray detector 
is capable of operating reliably on the Grand Tour spacecraft. 
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PART IV 

SUMMARY AND ACKNOWLEDGMENTS 



71 . 


The present best choice of instrument thus appears to be: 

1. Two separate multiwire proportional counters for 
redundancy. 

2. Passive collimation to restrict the field to ~5°, wire- 
grid modulation collimation to ~0. 1° angular resolution. 

3. No active shielding system around the counter body is 
required. 

4. Light passive shielding around any portion of the 
counter body exposed to space to absorb most of the 
cosmic x-ray background. 

Instruments of this description are well known in the laboratory 
and are widely used in x-ray astronomy on rockets and will soon be 
used On satellites such as OSO and HEAO. Development of the 
planetary x-ray experiment is therefore largely accomplished. 
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